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CsJ ■ ABSTRACT 

Studies of extinction curves provide insights into the properties of interstellar dust. Until recently, however, very few extinction curves 
'~, , existed outside the Local Group. GRB afterglows are well suited to extinction studies due to their brightness, simple power-law spectra 

• and their occurrence in distant star forming galaxies. In this paper we present results from the SED analysis of a sample of 41 GRB 

afterglows, from X-ray to NIR wavelengths. The sample is based on spectra from VLT-FORS, with additional data primarily from 
Swift. This is the largest sample of extinction curves outside the Local Group and, to date, the only extragalactic sample of absolute 
extinction curves based on spectroscopy. Estimates of the distribution of restframe visual extinctions, the extinction curves, and the 

O intrinsic spectral shapes of GRB afterglows are obtained. Their correlation with H I column density as well as total and gas-phase 
' metal column density are examined. The line-of-sight gas-to-dust and metals-to-dust ratios are determined and examined as a function 

\^ I of total column density, ISM metallicity and redshift. The intrinsic SEDs of the afterglows show that approximately half the sample 

require a cooling break between the optical and X-ray ranges. The broken power-law SEDs show an average change in the spectral 
index of A/? = 0.51 with a very small standard deviation of 0.02 (excluding the outlier GRB 080210). This is consistent with the 
expectations from a simple synchrotron model. Such a close convergence of values suggests that the X-ray afterglows of GRBs may 
^ ' be used with considerably more confidence to set the absolute flux level and intrinsic spectral indices in the optical and UV. Of the 

5-H , sample, 63% are well described by a featureless (SMC-type) extinction curve. Almost a quarter of our sample is consistent with no 

j/j ' significant extinction (typically Ay < 0.1). The 2175 A extinction bump is detected unequivocally in 7% of our sample (3 GRBs), 

' which all have Ay > 1.0, while one afterglow has a very unusual extinction curve with a sharp UV rise. However, we can only say 

that the bump is not present in about a quarter of our sample because of low extinction or lack of coverage of the 2200 A region. All 
the afterglows well fit with SMC type curves have moderate or low extinction, with Av < 0.65. This suggests that the SMC extinction 
curve is not as nearly-universal as previously believed and that extinction curves more similar to those found in the Galaxy and the 
LMC may be quite prevalent. We find an anti-correlation between gas-to-dust ratio and metallicity consistent with the Local Group 
relation; we find, however, no correlation between the metals-to-dust ratios and the metallicities, redshift and visual extinction; we find 
no strong conelation of the extinction column with metallicity either. Our metals-to-dust ratios derived from the soft X-ray absorption 
are always larger (3-30 times) than the Local Group value, which may mean that GRB hosts may be less efficient at turning their 
metals into dust. However, we find that gas, dust, and metal column densities are all likely to be influenced by photo-ionization and 
^SJ , dust destruction effects from the GRB to differing extents and caution must be used in extrapolating the ratios of dust and gas-derived 

. properties from GRB afterglows to the star-forming population in general. 
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1. Introduction optical and UV energy bands, revealing information about the 

. . ^, dust grain size and composition by consideiing the amount of 

^ Extinction of light is the scattering and absorption of photons ^. ^^ j^^^ scattering and absorption. The extinction curves 

■ ■ along the travel path from an object to the observer. Dust grams ^j^^ ^.j^ ^ ^^^^ Large and Small Magellanic 

are the major contributor to extinction in the optical and ul- cloud s (LMC and SMC) have been modeled in a variety o f 

traviolet (UV) range. More than 30% of the UV and opti- ^avs (iFitzpatrick & Massall986tlC^dem et al.ll989UPeirr99l . 

cal light from stars m the Universe may be absorbed and re- ^^^^^ extinction laws vary significantly, largely due to the pres- 
radiated thermally at infrared and millimeter wavelengths by the , , • > .c T 2 ■ ■ ^ , ■ , 

, • J , J7 — m Urvrvg L F T^"^ — r~n l-iXrvo k ence and relative strength of the 2 1 75 A extinction feature which 

obscuring dust (ILagache et al.M2005t IFranceschini et al.M2008D . ■ , . , • , ■• 

■ ■ — — T — j' — r~ ~~: — Z — . ■ — IT IS believed to be mainly due to absorption by graphite grains 

Extinction curves are the standard tool to study dust in the jrr — - — r— =; i m^^ LIH — ■ l-irvrvo k . 

(IStecher & Donnll9o5tlDramel2003D . Extinction curves are tvp- 

Send offprint requests to: tayyaba@dark-cosmology.dk ically obtained by comparing pairs of stars, however, outside the 

* Based on observations collected at the European Southern Local Group of galaxies, the method of estimating extinction 

Observatory (ESO) utilizing the 8.2m Very Large Telescope (VLT), curves from stellar atmosphere models is not applicable because 

Chile, under programs 075.D-0270 (PI: Fynbo), 077.D-0661 (PI: individual stars are difficult or impossible to observe. 
Vreeswijk), 077.D-0805 (PI: Taghaferri), 078.D-0416 (PI: Vreeswijk), 

079.D-0429 (PI: Vreeswijk), 080.D-0526 (PI: Vreeswijk), 081.A-0135 Long duration gamma-ray bursts (GRBs) are signposts of 

(PI: Greiner), and 281.D-5002 (PI: Delia Valle). star formation due to their association with the deaths of short- 
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lived, massive stars (e.g., WooslevI 19931: Galamaetal. 1998; 
Hioith et al.' '2003; 'Stanek et al.' '2003'; 'Delia Valle et al. '2003J 
Malesani et al. 2004; Campana et al. 2006; Starling et al. 201 1). 
GRBs are excellent probes to study extinction curves in dis- 
tant galaxies due to their intrinsic brightness, simple spectra and 
their occurrence in dense, star forming environments. GRB after- 
glows typically have intrinsically featureless piecewise power- 
law spectra from the X-ray to the optical/NIR and radio wave- 
lengths. Theoretically, the most promi nent explanation of GRB 



afterglows is the fireba l l m odel (e.g. iRees & Meszaros 



Meszaros & Reeslll997l: lOaia ma et all 1 1998b 
Piranlll999l: iGranot & Sarill2002t IZhang et al.l 



Sari et al 



1992 



1998 



20061) which pro- 



vides some notable successes in explaining the observations. 
According to the fireball model the GRB afterglow originates 
from synchrotron radiation produced by the interaction between 
the ultra-relativistic jet and the surrounding interstellar medium 
(ISM). Cooling of electrons produces a break in the synchrotron 
spectrum, which can sometimes be located between the X-ray 
and optical/NIR. The cooling frequency will lead to a well- 
defined change in spectral slope of A/3 - /St - /3\ - 0.5 resulting 
in a softer spectrum on the high energy side. Therefore, the X- 
ray-optical SED is expected to be a single or broken power-law. 
The simple spectral shape of GRBs and the lack of absorption 
above ~ 3 keV allow us to place constraints not only on the red- 
dening due to dust, but on the absolute extinction and thereby 
obtain absolute extinction curves for sources over a wide range 
of redshifts (z = 0.0086 - 8.2). 

In the era of the Swift satellite dGehrels et al.l |2004 and 
fast responding ground-based telescopes, better- sampled GRB 
afterglow lightcurves and SEDs enable us to get more com- 
plete extinction curves and even samples of extinction curves. 
To date, SMC-type extinction (i.e. without a 2175 A bump), 
has provide d the bes t repre s entation of GRB a fterglow ex- 
tinction (iBloometal.1 119981; I Jensen etal.1 120011; iFvnboetal 



200 It iPrice et al.ll200li ISavaglio & Fallll2004t Ijakobsson et al 



2004a iKann et alJl200a IStarling et al.ll2007l; iKann et alJl201(i 
Schadv etal.]|2010h . The 2175 A feature has been firmly de- 
tected in only two GRB afterglow s so far: GRB 070802 
(lElfasd ottir et al.l 
( Procha ska et al.l 



Kriihler et all 200: 



rr^-**r GRB 080607 

Perlev et alT 201 Ibb . Grey/wavelength- 
independent extinctio n laws have also been claimed for a 
few GRB a f terglo ws (iGalama et al.l [2003; Nardini et all 120081: 
iPerlev et"al] |2008|) . Dust with an unusual extinction curve, 
claimed to be due to super nova-origin dust, h as been suggested 
in two cases, GRB 071025 dPerlev et al.ll2010() and GRB 050904 
dStratta et al.ll2007h . but rul ed out for the latte r case with revised 
multi-epoch SED analysis dZafar et al.ll2010l) . 

In this paper, we present the SED analysis of a sam- 
ple of forty one GRB afterglows observed mainly with the 
Very Large Telescope (VLT). This is the largest study to date 
of extinction curves outside the Local Group of galaxies. It 
is also unique as the first spectroscopic sample of absolute 
extinction curves ( cf. K ann^et al. 2006; Starling et al. 2007^; 
Schadv et al.'2007b, 2010l: lKann et all2010l;lKruhler et air201U 



Gallerani et al.a2010 ). Our main objectives are: ;) to examine the 
extinction curves of star-forming galaxies at all redshifts, //) to 
understand the properties of the dust and their relationship to the 
gas and metals in these galaxies, ///) to infer information about 
the typical spectral shapes of GRB afterglows. 

To the best of our knowledge, this paper presents the first 
set of absolute extinction curves based on spectroscopy. Using 
spectroscopy in addition to photometry has significant benefits in 
such an analysis: /) the actual shape of the continuum and of the 



extinction curve can be determined, including the width, strength 
and central wavelength of the 2175 A bump where it is present 
;/) the effects of gas absorption on the spectrum can be removed, 
including Lycf at z > 2 and the Lya forest at higher redshifts, 
;//) precise redshifts can be obtained from the spectrum, iv) the 
potential influence of intervening absorbers can be determined. 

In §2 we describe the X-ray, UV/optical and NIR data re- 
duction. In §3 we describe the different models used to fit the 
data. The results of our modeling are presented in §4. A dis- 
cussion on the implications of our results is given in §5, and 
our conclusions are summarized in §6. A cos mology where 
Hp = 72kms'' Mpc ',Qa = 0.73 andfi,n = 0.27 dSpergel et al.1 
l2003h is assumed. Throughout spectral indices are denoted by y8 
where /3i is the low energy (typically the optical/NIR) slope and 
/?2 is the high energy (typically X-ray) slope. The temporal in- 
dices are denoted by a: Fy{v, t) oc r^v^^. All errors are Icr unless 
otherwise specified. For the cases where there is no detection, 2cr 
upper limits are reported. 

2. Multi-wavelength observations and data 
reduction 

2. 1 . Ensemble selection criteria 

The sample consists of the flux-calibrated spectra from 
iFvnbo etaP (|2009, hereafter F09) taken with the Focal Reducer 
and low dispersion Spectrographs (FORSl and FORS2) on 
the VLT with photometric observations in at least R or I- 
bands. Forty one GRB afterglows fulfilled our selection crite- 
ria (see Table [TJ. The spectra were obtained through target-of- 
opportunity (ToO) programs. 

2.2. Optical spectroscopy 

With the aim of constructing composite SEDs, we collected 
flux calibrated spectroscopic data for each GRB afterglow listed 
in Table [T] (see also F09). The sample includes GRB after- 
glows from April 2005 to September 2008. The longslit spectra 
were obtained with eithe r the FORS 1 or FORS2 spectrographs 
dAppenzeller et aT ][T998h using the 300V grism for all bursts ex- 
cept for GRB080319B (600B) and GRB 080913 (600z). This 
yields a spectral resolution of R - A/AA - 440, 780 and 1390 
in the 300V, 600B and 600z grisms, respectively, with a slit 
of 1". Further details of the optical spectroscopy are available 
in F09. The data were reduced using standard tasks in IRAF 
and corrected for co smic rays using the method developed by 
Ivan DokkumI d200ll) . All the spectra were flux calibrated using 
spectrophotometric standard stars observed on the same night of 
the afterglow spectra. All FORS 1/2 spectra from F09 are used 
except GRB 060526, GRB 060719, GRB 060807, GRB 07041 1, 
GRB 070508, GRB 070306, GRB 080411, GRB080413B and 
GRB 080523. There were no suitable standard star spec- 
tra obtained near the time of observations of GRB 060526, 
GRB 080411 and GRB 0804 13B. The spectra of GRB 060719, 
GRB 070411, GRB 070508 and GRB 080523 have extremely 
low S/N ratios. GRB 060807 has no redshift estimation. The 
spectrum of GRB 070306 was dominated by the host galaxy. The 
spectrum of GRB 080607 was obtained w ith Keck using th e Low 
Resolution Imaging Spectrometer (LRIS; IOke et al.ll995l) and is 
therefore formally outside our sample. We include it here as a 
comparison object, since it is one of only two known GRB after- 
glows with a confirmed 2175 A bump in the host galaxy prior to 
this sample. We used photometry from acquisition camera imag- 
ing, mostly performed using the R filter, to set the absolute flux 
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level of the spectra. The time of the acquisition camera image 
was then used as the mean time of the SED. The magnitudes 
of the afterglows are obtained using ESO zero-points from the 
night of the observation. 

To avoid contamination from absorption caused by the Lya 
forest, data (both spectroscopic and photometric) with restframe 
wavelength A < 1216 A were excluded from the SED fitting. 
Furthermore, emission and absorption lines arising from differ- 
ent metal species and atmospheric molecular absorption lines 
were removed from the optical spectra. 

2.3. Photometry 

The 30 cm Ultraviolet Optical Telescope (UVOT) onboard Swift, 
has six UV/optical filters: v, b, u, uvwl, uvm2 and uvw2 
jRoming et al.l 120051) . The data are made available at the UK 
Swift Science Data Centre. The UVOT filter lightcurves were 
extracted using the UVOT imaging data tool UVOTMAGHIST. 
The fluxes in each band were obtained by using the stan- 
dard UVOT pho tometric calibration techniques defined by 
iPoole et aP (l2008h . For the ground-based optical and NIR data, 
lightcurves for each filter were constructed using data either 
from refereed publications or from GRB Coordinates Network 
(GCN) circulars. For ease of access to the data published in GCN 
circulars for each GRB, the information was obtained via the 
GRBlo^ database ( lOuimbv et alj|2004l) . To generate SEDs at 
instantaneous epochs, the magnitudes in the relevant filters were 
obtained by interpolating or extrapolating the lightcurves to the 
time of interest. For this purpose, usually the /?-band lightcurve 
decay index is used for extrapolation and interpolation in other 
bands since it is typically the most complete. 

The optical spectra and optic al/NIR photom etry have been 
correc ted for Galactic dust extinction using the ISchlegel et alj 
119981) dust extinction maps. The E(B - V) values along the line 
of sight for each burst are given in Table [Tjand were retrieved 
from the NASA Extragalactic database web-calculator. 
Several wo rks have attempted to verifv the.Schlegel eFall (Il998h 
results (e.g. lDutra et al.l2003HSchlaflv et al I2010I) . finding over- 
all a good agreement at lo w values of A y . Som e uncertainty is 
still present, e.g. recently ISchlaflv et al.l (1201 Ol) suggested that 
these values might be systematically overestimated by ~ 14% 
on the average. Fo r our analysis we co rrected the Galactic 
E{B - V) values of ISchlegel etal] (Il998h for flie 14% overes- 
timation (see Table [T]l. However, extinction at high redshift has 
a much greater effect for a given dust column on the observed 
spectrum since we are observing dust effects in the restframe 
UV of the GRB host. For example, for a GRB host at z = 2 with 
SMC extinction, the effective observer frame is ~ 4 times 
larger than the equivalent Galactic Ar. Discrepancies between 
thelSchlafly et al. (2010) and Schlegel et al. ( 1998) values in the 
Galactic extinction are therefore small compared to the typical 
uncertainty on the measured restframe absorption and extinction 
values. 

2.4. X-ray data 

The GRB afterglows present in our sample hav e been all de- 
tected by Swift. The X-Ray Telescope (XRT; iBurrows et al.l 
l2005l) onboard Swift performed the X-ray observations in each 
case and the X-ray data for each GRB have been obtained 



' |http://grblog.org/grblog.php 
^ fhttp://nedwww.ipac.caltech.ei 



Table 1. Basic data on the flux calibrated spectroscopic sample. 
The details of the sample are provided in the columns as (1) burst 
name, (2) Galactic Hi column density, (3) time of SED, time 
since G RB trigger, (4) Galactic E{B - V) from the lSchlegel et al.l 
(Il998l) extinction maps corrected for 14% overestimation, (5) 
redshift, and (6) magnitude of the afterglow from the acquisition 
image. 
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from the UK Swift science data centefl The X-ray spectrum 
for each GRB afterglow was reduced in the 0.3-10.0 keV en- 
ergy range using the HEAsoft software (version 6.10). For each 
after glow, the X-ray lightcurv es were produced from the XRT 
data (lEvans et al. 1 120071 12OO9I) . To g et a better estimate of the 
X-ray flux, we fitted a decay model dBeuermann et al.ll 19991) to 
the X-ray lightcurves. For X-ray spectral analysis we used pho- 
ton counting (PC) mode data and source counts were extracted 
within a circular region centered on the source, and background 
counts from a circular region ~ 10 times greater than the source 

' Ihttp : //www . swift . ac . uk/ swi£t_portal/ archive . phpl 
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region. We used XSELECT (v2.4) to extract spectral files from 
the event data in the 0.3-10 keV energy band. To avoid pile-up 
all X-ray spectra were taken from time intervals where the mea- 
sured count rate is less than 0.6 counts s The X-ray spectral 
files were grouped to 12-20 counts per energy channel. Effective 
area files for each spectrum were created using the FTOOLS 
XRTMKARF recipe and bad columns were corrected by using the 
exposure maps. Response matrices (vll) from the Swift XRT 
calibration files were used. The X-ray spectra were obtained over 
an interval as close as possible to the mid-time of the optical 
spectra. The flux of the final X-ray spectrum was corrected to the 
SED time by taking the ratio of the X-ray lightcurve model fits 
at the photon weighted mean time (PWMT) of the X-ray spec- 
trum and the SED time . The spectra were fitted within XSPEC 
(vl2.6.0: lArnaudl[T996h with a Galactic-absorbed power-law, 
with absorptio n from the Galactic ne utral hydrogen column den- 
sity taken from lKalberla et al.l(l2005l using the nH FTOOL), and 
additional absorption from the GRB host galaxy. The Galactic 
column density value used for each GRB is given in Table|5] The 
equivalent neutral hydrogen column density from the host galaxy 
of each GRB was estimated from the soft X-ray absorption and 
is denoted here as A^h.x- A^h.x is m odeled with X SPEC where so- 
lar abundances were assumed from lAsplund et al.i 620061) . A^h.x, 
presented in units of equivalent hydrogen column density for a 
solar abundance of the elements, is essentially a direct measure 
of the metal column density. 



3. Data analysis 

We describe the afterglow continuum emission with single or 
broken power-law models. For each burst, we tried both laws to 
fit the observed data, where the spectral break must be in the 
observing window (10'"^Hz< Vbieak ^ 10'^ Hz). On top of the 
power-law fits, we added empirical extinction laws correspond- 
ing to die SMC iP ei 1992), MW jFitz patrick & Massa 198^ 
ICardelU et al.1IT98l . or LMC (lPeilll992l) (see ED for details), 
in order to account for the effect to dust in the GRB local envi- 
ronment. The X-ray to optical/NIR GRB afterglow SEDs were 
generated at specific epochs (the time of the acquisition image 
of the spectrum), and modeled in IDL. To get errors on each 
parameter we used 1000 Monte Carlo (MC) realizations of the 
data, with the mean value of each datapoint set to the observed 
value and a gaussian distribution with a width corresponding to 
the Icr measured error on that datapoint. For each simulated set 
we fitted the function and estimated the best fit parameters. The 
error for each parameter was calculated using the standard devi- 
ation of this distribution. To determine whether the use of a bro- 
ken power-law was required over a single power-law, we used 
an F-test probability < 5% as the threshold. The results from the 
spectral analysis of our spectroscopic GRB sample are provided 
in Table|2] The number of degrees of freedom in Table|2]is calcu- 
lated by taking the number of pixels in the optical spectra. Most 
of the spectra are taken with FORSl/2 grism 300V yielding four 
pixels per resolution element. 

3.1. Dust models 

The dust extinction was modeled with empirical extinctio n laws 
dFitzpatrick & Massalll986l:ICardeni et al.ll 9891: [PeiFl 9921) to de- 
termine the absolute extinction and the characteristics of the 
extinction. From the ,Pe i ( 199Z) extinction laws we used only 
the SMC extinction law. He reafter, we will refer to the se three 
extinction models as FM dFitzpatrick & Massal Il986h . CCM 



Table 2. Results of fits to the SEDs. The columns give the burst 
name, and then for each model the xl, the number of degrees of 
freedom (dof) and the null hypothesis probability (NHP) for the 
fit are provided. The F-test probability is computed for each SED 
comparing the single (PL) and broken power-law (BPL) models. 
The best fit models are denoted f. 



SMC 

PL BPL 



GRB 


Xl (dof) 


NHP% 


Xl (dof) 


NHP% 


F-test prob. 


050401 


1.09 


(1011) 


2.0 


1.08 


(1009)+ 


2.6 


0.02 


050730 


1.23 


(659) 


0.003 


0.98 


(657)+ 


60 


< 0.01 


050824 


1.13 


(958) 


0.2 


0.97 


(956)+ 


71 


< 0.01 


060115 


1.15 


(1014)"'" 


0.06 


1.15 


(1012) 


0.05 


1.00 


060512 


1.02 


(1452) 


33 


0.93 


(1450)+ 


96 


< 0.01 


060614 


1.01 


(1254) 


37 


0.82 


(1252)+ 


100 


< 0.01 


060707 


0.98 


(1088)"f 


65 


0.98 


(1086) 


64 


1.00 


060708 


0.86 


(867) 


100 


0.85 


(865)+ 


100 


< 0.01 


060714 


1.02 


(1043) 


35 


0.97 


(1041)+ 


72 


<0.01 


060729 


0.99 


(1140)t 


54 


0.99 


(1138) 


53 


1.00 


060904B 


0.80 


(1198)+ 


100 


0.80 


(1196) 


100 


1.00 


060906 


1.01 


(786) 


42 


0.94 


(784)+ 


87 


< 0.01 


060926 


0.95 


(907)t 


86 


0.95 


(905) 


85 


1.00 


060927 


1.03 


(274)t 


35 


1.03 


(272) 


32 


1.00 


061007 


0.97 


(1341)+ 


78 


0.97 


(1339) 


77 


1.00 


061021 


1.26 


(1703) 


0.0 


0.98 


(1701)+ 


70 


<0.01 


061110A 


0.98 


(1093)+ 


63 


0.98 


(1091) 


61 


1.00 


061 HOB 


0.95 


(896) 


85 


0.94 


(894)+ 


88 


0.02 


070110 


1.09 


(1087) 


0.17 


1.04 


(1085)+ 


18 


< 0.01 


070125 


1.03 


(790) 


25 


0.97 


(788)+ 


72 


< 0.01 


070129 


0.93 


(1017) 


94 


0.91 


(1015)+ 


98 


< 0.01 


070506 


0.99 


(775) 


57 


0.98 


(773)+ 


65 


< 0.01 


070611 


0.77 


(1248)+ 


100 


0.77 


(1246) 


100 


0.77 


07072 IB 


1.00 


(830)+ 


47 


1.00 


(828) 


46 


0.76 


071020 


1.09 


(1020) 


2.1 


1.06 


(1018)+ 


10 


< 0.01 


071031 


2.13 


(614) 





1.05 


(612)+ 


17 


< 0.01 


071112C 


1.67 


(1287) 





1.00 


(1285)+ 


48 


< 0.01 


071117 


n net 
0.99 


11/1) 


56 


0.98 


,1 lo9) 


71 


< 0.01 




1.10 


(964) 


1.7 


1.00 


(962)+ 






0803 19B 


0.83 


(1468)+ 


100 


0.83 


(1466) 


100 


1.00 


080520 


1.00 


(1028)+ 


54 


1.00 


(1026) 


51 


1.00 


080707 


1.07 


(743) 


9.3 


1.03 


(741)+ 


30 


<0.01 


080721 


0.96 


(731)+ 


76 


0.96 


(729) 


74 


1.00 


080905B 


0.98 


(1029)+ 


68 


0.98 


(1027) 


66 


1.00 


080913 


0.58 


(81)+ 


100 


0.58 


(79) 


100 


1.00 


0809 16A 


0.93 


(1170) 


96 


0.91 


(1168)+ 


99 


< 0.01 


080928 


0.99 


(1157)+ 


57 


0.99 


(1155) 


55 


1.00 








FM 












PL 






BPL 






070318 


1.15 


(1392) 


0.001 


0.88 


(1390)" 


100 


< 0.01 


070802 


0.54 


(1381)+ 


100 


0.54 


(1379) 


100 


1.00 


080605 


0.86 


(859)+ 


100 


0.86 


(857) 


100 


1.00 


080805 


0.93 


(1391)+ 


98 


0.93 


(1389) 


97 


1.00 


080607 " 


1.79 (39)+ 


0.19 


1.79 (37) 


0.17 


1.00 



" Keck spectrum. 



dCardelli et al.ll 19891) and SMC (lPeilll992l) . The important char- 
acteristic of the SMC law used here is the lack of the graphite 
2175 A b ump and the even steeper rise in the far UV, rising faster 
flian 1 /A dPrevot et al.lll984t iGordon et al.ll2003l) . ICardeUi et al.] 
dl989h showed that Galactic sightlines could typically be well- 
fit with a parameterization of the extinction curve that relied 
only on a single parameter, the ratio of absolute to selective ex- 
tinction, Rv - Av/E{B - V). Comparable to CCM, FM pro- 
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vides greater freedom in reproducing the extinction curves . To 
investigate restframe visual extinction, we fit all three extinc- 
tion laws for each case and found that the best fitting results and 
null hypothesis probabilities for SMC and FM were similar in 
most cases. We therefore selected the SMC as by far the simpler 
model and used it to characterize the extinction in 37 cases (see 
Table |2]i where the null hypothesis probability for the fit with 
FM was not significantly better. The FM parameterization was 
used for the cases where the SMC could not provide a good fit 
because of the presence of the 2175 A bump (see Table|2]i. The 
CCM curve yielded a poorer fit to the observed data in every 
case (see below) 

The continuum emission from an afterglow is bel ieved to be 
dominated by synchrotron radiation (ISari et al.lll998l) described 
by a power-law given as 



Fqv 



(1) 



where Fq is the flux normalization, v is the frequency and /? is 
the intrinsic spectral slope. The observed flux is extinguished 
due to dust in the Milky Way, and in objects between our Galaxy 
and the GRB, the latter typically dominated by the G RB host 
galaxy (though see GRB 060418: lEUison et ail l2006l) . As ex- 
plained above (Sect. I2.3l l, we correct ed the optical/NIR d ata for 
Galactic foreground extinction maps (ISchlegel et al.ll998h . With 
knowledge of the unextinguished flux from X-rays and the red- 
shift, the extinction corrected spectral slope and restframe Ay 
can be obtained. Therefore the extinguished flux will be given as 



-0.4A., 



(2) 



Here is the extinction in the host galaxy of the GRB as a func- 
tion of wavelength A. The wavelength dependence on dust ex- 
tinction observed in SMC, LMC and MW environments is well 
reproduced by the standard dust models (iFitzpatrick & Massal 
[l986; .Cardelli et al. 1989; Pel 1992 .) For S MC, LMC and MW 
the extinction law introduced by |Pei d 19921) is given as 



(A/Ai)"' + (Ai/A)"^ + bi 



(3) 



For SED modeling we used the SMC parameter ization where 
Rv = 2.93 for the SMC extinction law (|Peilll992h . The param- 
eters a,;, bj and «, are different for each extinction model (see 
lPeilll992l for the parameter list). A param eterization for extinc- 
tion in the Milky Way was advocated bv ICardelli et al.l (Il989h . 
It provides a better estimation of the far UV curvature with fixed 
bump strength and varying Ry- It is defined as 



Av\a{x) + 



b(x) 
Rv 



where x - A 



(4) 



a(x) and b(x) are wavelength dependent polynomials controlling 
the UV, optical and NIR regimes of the extinction curve. Finally, 
we fitted the data with an FM extinction model. The parame- 
terization for the FM extinction curve in the UV range is given 
as 

A,^Avl-^kiA-V) + l] (5) 



where 



k(A - V) 



Rv 



c\ + cox + c^Dix, X(), y) 

ci + C2X + ct,D(x, xo, y) + C/i{x — c^Y 



X < C5 
X > C5 



and the Lorentzian-like Drude profile, which represents the 
2175 A bump, is given as 

D(x, xo, y) - 



(x^ - Xq)2 + x^y^ 

The underlying UV linear component is specified by the in- 
tercept, ci, and the slope, C2. The ~ 2175 A absorption bump 
is controlled by C3 (bump strength), xq (central wave number), 
and y (the width of the bump). The parameters q and C5 give 
the far-UV curvature. The extinction properties in the NIR and 
optical are determined us ing spline interpolation points (see 
IFitzpatrick & Mas^l2007h . The FM parameterization provides 
greater freedom for fitting the extinction curve; in particular, it 
allows characterization of the curve, for example in the strength, 
width and position of the 2175 A bump. 



4. Results 

4. 7. The SED of the afterglow 

GRB afterglows are well described by a coherent synchrotron 
emission from accelerated electrons with a power-law distribu- 
tion of energies promptly behind the shock. These electrons then 
cool both adiabatically and by emitting synchrotron and inverse 
Compton radiation. At late times, the electron population accel- 
erated by shock me chanisms is exp ected to radiate in the slow 
cooling regime (see lSari et al.lll998b . The GRB spectrum in that 
regime is defined by a broken power-law model. The transfer to 
the slow cooling regime causes a cooling break in the intrinsic 
spectrum, typically occurring between the X-ray and the opti- 
cal/NIR wavelengths. The broken power-law model is given as 



" hrenk 



V < Vbreak 

V > Vbreak 



(6) 



where /3i and (^2 are the slopes of the optical and X-ray seg- 
ments in our model respectively. The parameter Vbreak is the 
cooling break frequency joining the two power-law segm ents. 
In accordance with the theory outlined bv lSari et al.l (1 19981) . this 
additional break in the spectrum will give a spectral change of 
/?! — P2 - AyS, where AyS — 0.5. In our SED modeling we treat 
Ay6 as a free parameter by allowing both slopes to vary indepen- 
dently in order to test the compatibility of the synchrotron model 
with our observed SED sample. From our sample, 21 SEDs have 
a spectral break between the X-ray and optical wavelengths. We 
find that 95% of the spectral breaks have a spectral change of 
Ay6 ~ 0.5, in agreement with the synchrotron emission model. 
The mean A/S of all 21 GRBs is 0.54 with a standard devia- 
tion of 0.13. This large standard deviation is dominated by a 
single outlier, GRB 080210, which has a A^S = 1.1 + 0.1. Of 
the 21 GRBs, GRB 080210 is the only case that deviates from 
Ay8 = 0.5 (see iDe Cia et al.ll201ll for a detailed analysis of the 
SED of GRB 080210). Considering only the remaining SEDs, 
(Ay6) = 0.51 with a st andard deviation of 0.02. This agreement 
with the prediction of ISari et"al ] (119981) is remarkable, and indi- 
cates that for the vast majority of SEDs taken at similar epochs, 
it is valid to assume that the amplitude of the spectral break is 
equal to exactly 0.5. In the future, for broad-band SED analysis, 
this should remove a degree of freedom in deriving extinction 
curves using GRBs, making the process more robust. 

This result provides strong support that the afterglow emis- 
sion is indeed synchrotron radiation, and that both the optical 
and X-ray emission are produced by the same mechanism in the 
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Fig. 1. A/3 against the X-ray spectral slope. The dashed line cor- 
responds to = 0.5. GRB 080210 is an outlier and has a softer 
X-ray slope compared to the rest of the sample. 

vast majority of cases. While this is a key element of the fireball 
external shock model, our findings are not in complete agree- 
ment with it: we checked whether our data satisfy the closure re- 
lations be t ween the spectral and temporal indices as defined by 
ISari et al l (ll998h . In most cases, we found that the observed tem- 
poral indices were not consistent with those expected from the 
fireball model, even though the specific prediction of Afi = 0.5 
does hold. This is not surprising. The synchrotron emission pro- 
cess is only one of the elements of the fireball model. The clo- 
sure relations depend not only on the spectral shape, but also on 
the temporal evolution of the fireball properties (radius, Lorentz 
factor, magnetic field), on the assumption that the microphysics 
parameters and sb do not vary in time, and on several other 
details. The shape of the synchrotron spectrum, on the contrary, 
is not sensitive to these additional parameters, and therefore the 
prediction that A/3 = 0.5 is much more robust. Indeed, our analy- 
sis is not dependent on the fireball dynamics, but it just assumes 
(and confirms) that the optical and X-ray emission are due to 
synchrotron emission. 

4.2. Dust extinction 

The restframe visual extinction. Ay, is determined for the forty 
one afterglows in our sample, plus GRB 080607. Approximately 
two thirds (26/41), have low/moderate extinction values with 
Ay ranging from 0.06 ± 0.02 to 0.65 + 0.04, and are well- 
fit with SMC-type extinction. The mean visual extinction for 
these 26 GRB afterglows is (Ay) - 0.24 with a standard de- 
viation of 0.14. Roughly one quarter of the sample (11/41) is 
consistent with no dust reddening, with a mean 2cr upper limit 
(Ay) < 0.10 assuming an SMC type dust extinction curve. More 
than 7% of the full sample (~ 9% of the extinguished after- 
glows) are found to have a well-defined bump consistent with the 
2175 A feature. Of these bursts, GRB 070802 was known before 
dlllasdottir et al. 2009; Kruhler et al. 2008), while GRB 080605 
and GRB 080805 are new (but see F09), doubling the number 
of such bumps unequivocally known in GRB host galaxies. We 
also include the other definite case of a bump in a GRB host 
extinction curve, GRB 080607, in our analysis as a comparison 
extinction curve. All of these cases have Ay > 1, and are by far 
the most extinguished bursts in our sample, with a mean extinc- 
tion more than five times the average extinction of extinguished 



Table 3. Param eters of the best fit with the SMC extinction curve 
(Rv = 2.93) of |Pe3 (Il992h . The values of /?2 and Vb,eak are pro- 
vided for the cases where a broken power-law is the best fit. 



GRB 


/3i 








82 






log Vb 

Hz 






Av 
mag 


050401 


0.39 ± 


0.05 





89 


± 





08 


16.32 + 





89 


0.65 + 0.04 


050730 


0.16± 


0.02 





66 


± 





02 


16.00 ± 





17 


0.12 + 0.02 


050824 


0.40 ± 


0.04 





92 


± 





07 


15.20 + 





33 


0.15 + 0.03 


060115 


0.77 ± 


0.03 


















0.10 + 0.02 


060512 


0.68 ± 


0.05 


1 


16 


± 





10 


15.37 + 





45 


< 0.12 


060614 


0.47 ± 


0.04 


1 


00 


± 





05 


16.06 + 





21 


0.11+0.03 


060707 


0.59 + 


0.02 


















0.08 + 0.02 


060708 


0.67 ± 


0.04 


1 


16 


± 





11 


16.92 + 


1 


04 


0.14 + 0.02 


060714 


0.44 ± 


0.04 





98 


± 





11 


15.86 + 





53 


0.21 + 0.02 


060729 


0.78 + 


0.03 


















0.07 + 0.02 


060904B 


0.90 ± 


0.02 


















0.34 + 0.03 


060906 


0.56 ± 


0.02 


1 


10 


± 





10 


17.10 + 





14 


< 0.09 


060926 


0.82 ± 


0.01 


















0.32 + 0.02 


060927 


0.86 ± 


0.03 


















< 0.12 


061007 


0.78 ± 


0.02 


















0.34 + 0.03 


061021 


0.55 ± 


0.02 


1 


07 


± 





04 


16.93 + 





24 


< 0.10 


061110A 


0.79 ± 


0.02 


















< 0.10 


061 HOB 


0.73 ± 


0.03 


1 


23 


± 





11 


16.90 + 





15 


0.23 + 0.03 


070110 


0.66 ± 


0.02 


1 


15 


± 





10 


16.84 + 





18 


< 0.10 


070125 


0.55 ± 


0.04 


1 


06 


± 





07 


14.95 + 





48 


0.30 + 0.04 


070129 


0.58 ± 


0.02 


1 


12 


± 





09 


16.98 + 





13 


0.28 + 0.02 


070506 


0.51 ± 


0.06 


1 


02 


± 





08 


16.00 + 





69 


0.44 + 0.05 


070611 


0.95 ± 


0.02 


















0.06 + 0.02 


07072 IB 


0.88 ± 


0.02 


















0.20 + 0.02 


071020 


0.15 ± 


0.05 





69 


± 





03 


14.79 + 





25 


0.43 + 0.04 


071031 


0.64 ± 


0.01 


1 


13 


± 





11 


15.35 + 





55 


< 0.07 


071112C 


0.37 ± 


0.02 





86 


± 





06 


15.67 + 





22 


< 0.08 


071117 


0.41 ± 


0.03 





94 


± 





08 


16.95 + 





92 


0.28 + 0.02 


080210 


0.25 ± 


0.03 


1 


36 


± 





14 


16.03 + 





15 


0.33 + 0.03 


0803 19B 


0.77 ± 


0.02 


















< 0.11 


080520 


0.78 ± 


0.02 


















0.22 + 0.02 


080707 


0.71 ± 


0.03 


1 


22 


± 





05 


16.90 + 





16 


< 0.12 


080721 


0.68 ± 


0.02 


















< 0.12 


080905B 


0.80 ± 


0.07 


















0.42 + 0.03 


080913 


0.79 ± 


0.03 


















0.12 + 0.03 


0809 16A 


0.49 ± 


0.05 





99 


± 





09 


17.00 + 





88 


0.15 + 0.04 


080928 


1.08 ± 


0.02 


















0.29 + 0.03 



bursts in the sample. Finally, GRB 070318 was not well fit by 
the SMC curve and does not have a 2175 A bump. Its extinction 
curve seems very unusual and is the subject of a separate investi- 
gation (Watson et al. in prep.). The Maiolino et al. (2004) extinc- 
tion curve associated with dust from supemovae was not used, 
but we infer from the good fit obtained with the SMC extinction 
curve in all cases without a 2 175 A bump, that this extinction 
curve is not required to fit the data. 

5. Discussion 

5.1. Tlie2175A feature 

The 2175 A absorption feature, as mentioned above, has been 
robustly detected in the optical spectra of four GRB afterglows 
(see Table |4]i. The feature is a broad dip centere d at 2175 A du e 
to excess extinction and was first discovered bv lSteched(ll965h . 
Notable characteristics are the width (~ 1 /im ') and its vari- 
ability (~ 20%), contrasted to the tight observed central wave- 
length of the bump (to within ~ 10 A, i.e. < 1%). Several candi- 
dates have been considered to explain the 2175 A b ump, ranging 
from carbonaceous materials (iHenrard et al.lll997l) to iron poor 
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Table 4. Best fit extinction curve parameters for the afterglow cases where FM is the best model (see Table. |2]i. GRB 080607 was 
observed with Keck telescope and is not a part of our spectroscopic sample. 



Parameter 


070802 


080605 


080805 


080607 


C\ 


0.02 + 0.25 


-0.81 ±0.16 


0.26 ±0.15 


1.19 ±0.43 


C2 (jum) 


0.82 ±0.16 


0.86 ± 0.07 


0.57 ± 0.06 


0.37 ±0.10 




1.46 + 0.29 


0.35 ± 0.04 


0.89 ±0.11 


2.96 ± 0.66 




0.26 ± 0.08 


0.26 ±0.10 


0.15 ±0.14 


0.17 ±0.10 


C5 (pm"') 


5.03 ±0.31 


6.30 ±0.41 


6.50 ±0.14 


5.88 ±0.34 


J (yum"') 


0.99 ±0.01 


0.90 ± 0.02 


0.90 ± 0.02 


1.20 ±0.04 


Rv 


T Q 1 +0.67 


9 Q^-0.17 


2.48 ± 0.39 






4.64 ±0.01 


4.55 ±0.01 


4.57 ± 0.01 


4.52 ± 0.02 


Pi 


0.88 ± 0.06 


0.76 ± 0.03 


0.81 ± 0.02 


96+" "' 


Av (mag) 


1.19 ±0.15 


1 20+""'' 


1.53 ±0.13 





12 



10 



G 
3 
z 



Measured Ay + 2a Upper limits 
Measured values 



n 



n 



n 



tinction curves toward the SMC, and one in fact does have a 



0.0 



0.5 



1.0 1.5 
Av (mag) 



2.0 



2.5 



Fig. 2. Distribution of the restframe visual extinction. Ay, of 
our flux calibrated spectroscopic sample. The solid curve cor- 
responds to detections; the dot-dashed curve to detections and 
2cr upper limits. 



silicate grains in th e form of partially hydrogenated amorphous 
Mg2Si04 particles dSteel & D uley 19871). The dominant hypoth- 
esis is some form of carbonaceous material in small grains (< 
50 A) , with polycyclic a romatic hydrocarbons ("Dulev & Seahr; 



1998t lDulevll2006t ICe cchi-Pestellini et al. 2008; Duley 200 



( Stecher & Donni,1965L .Mathis et al. 


1 19771 Draine & Led 


1984 


Drainelll989l ISorrelllll990l iMathisI 


1994 


Roulea 


u et al.l 


1997 


Will & Aannestad' '1999^, 'Andersen et al. 


2002; 


Clavton et al 



20031: iFitzpatrick & Massa 2007j). The absence or weakness of 



the bump in SMC and LMC sightlines can be explained by a dif- 
ference in the relative abundances of graphite and silicate grains 
(|Peilll992l:|Prevot et al.lll984 . The presence of the 2175 A bump 
in GRB afterglow SEDs indicates that these GRBs occur in more 
evolved and massive galaxies and have dust properties similar to 
our local disk galaxies. Beyond our Local Group these distant 
GRBs can provide information on the graphite content of small 
grains in the surrounding circumburst environments. 

It is notable that less than 9% of the GRB extinction curves 
show an unequivocal 2175 A bump. This is a small fraction, 
consistent with the prevailing notion that strongly star-forming 
galaxies, and in particular GRB hosts, have featurele ss dust 



19941 Schadv et al. 


l2010HKann et alJ 


2010l;rKruhler et alj|2011l; 


Greineret al.i i2011 


). However, there are only five known ex- 



2175 A bump ( Peil 19921) . Furthermore, the starburst sample of 
ICalzetti et al.l d 19941) with SMC-like attenuation curves is very 
small, containing only six objects. Therefore we examined our 
sample to see in how many SEDs a bump could have been de- 
tected. We found that for 12 cases the bump region is not signif- 
icantly covered by our spectra or photometry and in a few more 
the signal-to-noise ratio is very low. This is approximately 40- 
50% of the extinguished sample. It could be argued therefore, 
that in only about a quarter of the sample are we reasonably cer- 
tain that no bump is present. Furthermore, there is very clear 
evidence in this sample that afterglows with a 2175 A bump are 
preferentially heavily extinguished (Ay > 1). 

It is apparent that our spectroscopic sample is heavily reliant 
on optically bright afterglows; the bulk of bursts in the sample 
have an average magnitude of 2 1 . 1 with a standard deviation of 
1.7 magnitudes. But a MW style extinction gives 2-3 magni- 
tudes of extinction in the BVRI bands for Ay - 1 at z = 2. A 
restframe Ay > 1 would make spectroscopy unfeasible for the 
overwhelmingly majority of these bursts and a redshift would 
almost certainly not have been found for any except the bright- 
est few (apart from emission lines from the host galaxy). It is 
obvious that we are heavily biased against these high extinction 
sightlines. F09 report a dark burst fraction of 25-42% for that 
sample. If most of those dark bursts are similar to the most heav- 
ily extinguished bursts in our sample, many, if not all, may have 
a 2175 A bump in their extinction curves. Our conclusion is that 
the apparent prevalence of the SMC extinction curve in GRB 
hosts is very likely an artifact of the bias in favour of low ex- 
tinction sightlines and the frequent lack of good coverage of the 
bump region. We can only say that the bump is not present in 
about a quarter of our sample, and that it is present in 7%, and 
that we know there is a clear preference for the bump to be de- 
tected in the most extinguished sightlines. Currently it is there- 
fore not possible to say whether a featureless extinction curve 
or one with a bump is more common in GRB host galaxies. But 
it seems likely that the 2175 A bump is far more prevalent than 
previously believed. 



5.2. Constraints on Ry 

We attempted to constrain Ry specifically for the extinction 
curves that were well fit with SMC-type extinction. We fitted 
the data again with the FM model, this time fixing the FM pa- 
rameters to the SMC average extinction curve values given in 
iGordon et al. (2003) and fixing A/3 to 0.5 for the 20 relevant 
spectral break cases (see 34. Il l, but allowing Ry to vary freely. 
In most cases Ry is degenerate with Ay, so that we cannot con- 
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Table 5. Table listing the spectroscopic GRB sample with their 
name, Galactic column density, neutral hydrogen column den- 
sity when available from the literature, equivalent column den- 
sity (A^H.x) determined from the X-ray spectral fits and the 
Nh,x/Av ratio. We reported 2cr upper limits for A^h.x where de- 
tection is < 2cr significant. 
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strain Ry in any useful sense. For GRB 071020 and 080928 we 
can constrain Ry to 3.0^yg and 2.0;^q^ respectively (Icr errors 
for two parameters of interest). From these cases, we can deter- 
mine that it is the strong constraint on the absolute optical flux 
level that allows the constraint on Ry- In the former case it is 
because the spectral break is very close to the optical, and in the 
latter because the X-ray and optical fluxes are so similar In both 
cases the extinction is relatively high - Ay = 0.43 + 0.04 and 
0.29 + 0.03 respectively - and the X-ray spectrum fairly well- 
defined. To fix Rv in future observations will require that the 
unextinguished power-law is well-constrained, either through in- 
frared spectroscopy or better X-ray observations. The ideal way 
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Fig. 3. Ay against spectroscopically measured Hi column den- 
sity in log-log space. A^hi values are obtained from the literature 
and given in Table|5] The dashed line corresponds to the Nui/Ay 
ratio f or the SMC environment given by iWeingartner & Drairiel 
(I2001h . 



to do this is to use mid-infrared data in concert with good X-ray 
data. To date this has been published for only one GRB after- 
glow using observations with Spitzer (Heng et al. 2008), albeit 
with a low S/N result. While more mid-infrared observations are 
now unlikely, and better X-ray observations will not be routinely 
available in the foreseeable future, near-infrared spectra are now 
regularly obtained with the X-shooter spectrograph on the VLT 
and may allow us to obtain reasonably well-constrained Ry val- 
ues for GRB extinction sightlines in the next few years. 

5.3. Gas-to-dust ratios 

The H I column densities for these lines of sight have been mea- 
sured from the 1216 A Lya absorption line in F09. The ratio of 
A^Hi to Ay gives, in principle, an estimate of the gas-to-dust ratio 
in the host galaxy of the GRB (Fig. [3). However, the potential 
ionizing and destructive effect of the GRB and its afterglow may 
substantially complicate this picture. The vast majority of data 
points lie at values greater than the SMC gas-to-dust ratio - the 
dashed line in Fig. [3] corresponds to the gas-to-dust ratio of the 
SMC er ivironment with a valiie of N m/Ay - 1.6 x 10^^ cm"- 
mag ' (IWeingartner & Drainell200TI) . The only clear outlier is 
GRB 071021, with a very small neutral hydrogen column den- 
sity quoted in F09. However, the burst is at z = 2.1, and the 
spectrum is cut off in the blue. On close re-inspection of the 2D 
spectrum, it seems unlikely that such a strong limit can indeed 
be placed on the Lya line in the afterglow of this burst. Indeed, 
it seems possible that this burst could have a large enough neu- 
tral column density in hydrogen to make it perfectly consistent 
with the SMC value. Therefore, the only burst in our sample 
that securely yields a gas-to-dust ratio lower than the SMC is 
GRB 070802, with a gas-to-dust r atio consistent w i th Ga lactic 
sightline values. As pointed out bv lElfasdottir et akl (|2009|) . this 
afterglow sightline is highly unusual among GRBs in this re- 
gard, a fact strengthened by our results. However, it is clearly 
not a property of all sightlines with the 2175 A bump, since 
GRB 080607 has a gas-to-dust ratio a factor of ten higher Due 
to the low redshift, we do not have Hi column density mea- 
surements for GRB 080605 and GRB 080805. The overafl dis- 
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log N„, + [M/H] 



Fig. 4. Visual extinction against U V metal abundances in solar 
units. The metallicities obtained from the literature based on Zn 
element are marked with red stars. The rest of the data points are 
based on Si element metallicities (see Table |6l). GRB 080607 is 
indicated with a red square. The metallicity or Ay upper limits 
are represented by small circles or stars. The dashed line indi- 
cates metals-to-dust relation for the Local Group (indicated as 
LG). 



tribution suggests that gas-to-dust ratios in GRB environments 
are typically SMC-like or higher. An obvious idea is that metal- 
licities that are SMC-like or lower are responsible for the high 
gas-to-dust ratios. However two arguments point against such a 
conclusion: first, the metals-to-dust ratio, where the metal col- 
umn density is derived from soft X-ray absorption is typically 
an order of magnitude or more above the values from the Loca l 
Group (see below and IWatson et al.l 120071: ISchadv et al.ll2010h : 
second, when we correct the H i column for the metallicity in 
Fig. m we show that, even corrected for metallicity, most GRB 
sightlines have a super-Local-Group gas-to-dust ratio. 

In Fig. ID we check if metallicity effects alone could be re- 
sponsible for the high gas-to-dust ratios without reference to the 
X-ray absorption. We derived a metal column density in solar 
units using log (A^Hi/cm"^) + [M/H] (equivalent to A^h.x)- These 
are metal column densities in the gas phase, derived from low 
ionization species, based on Zn and Si. The issue of depletion 
of the elements into dust grains is important. Zn and S are non- 
refractory elements and are the best tracers of the metallicity. We 
have used Zn metallicities wherever possible, however in many 
cases we have used Si to measure metallicity. Si is affected by 
depletion of the metal out of the gas phase, but less than strongly 
refractory elements like Fe an d Cr (e.g., iLedoux et al.l 120021 
iDrainelUOOSllFvnbo et al.ll2010l) . The [M/H] values are given in 
Tableland the A^hi values for the corresponding GRBs are taken 
from Table[5] In Fig.|4]we see that the low ionization metals-to- 
dust ratios are closer to the Local Group value than the gas-to- 
dust ratios, but are still typically higher. Conclusions on the gas- 
to-dust ratio of the host galaxies of GRBs are inevitably compli- 
cated by the fact that we know a majority of the gas along the 
line of sight to a GRB is i onized close to the burst (IWatson et al 



2007|:ISchadv et al.l20I0l:[Prochaska et al.l2007HVreeswiik et al 



2007; Ledoux et al ] |2009l) . and it is suspected that any dust near 
the GRB will be either sublimated or cracked by the GRB (e.g. 
^ruchter et al. 2001; Perna & Lazzati 2002). The X-ray absorp- 
tion, by contrast, is far less affected by photo-ionization, and so 
GRBs show much larger metals-to-dust ratios (Fig.|5). The most 
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Fig. 5. Ay measurements versus A^h.x for the spectroscopic GRB 
sample determined from X-ray spectral analysis and SED fit- 
ting. The dotted curve represents the metals-to-dust ratio for 
the Local Group environments (indicated as LG) assuming so- 
lar abundances. 



probable interpretation is that the column has been strongly af- 
fected by the GRB/afterglow, and the ionization distance and the 
dust sublimation distance are different. Therefore the metals-to- 
dust and gas-to-dust ratios derived from the optical/UV are not 
representative of the intrinsic ratios for the GRB host, but instead 
are due to the interaction of the distribution of the gas and dust 
along the column and the luminosity and spectrum of the GRB 
and its afterglow. Thus, the SMC-like and higher gas-to-dust ra- 
tios found in GRB afterglows, even after accounting for metal- 
licity, seem to be due to the competing effects of gas ionization 
and dust destruction. The ratios found are closer to the values in 
the Local Group than are the metals-to-dust ratios derived from 
X-rays. An obvious explanation for this is that since we know 
that hydrogen ionization certainly occurs in the surroundings of 
the GRB, the H i column densities are simply lower limits to the 
actual gas column. So far we have no unequivocal evidence for 
dust destruction in GRB surroundings, but it seems probable that 
it occurs; what we can at least say from these results is that hy- 
drogen ionization in the GRB surroundings is a more efficient 
process than dust-destruction, drawing the gas values closer to 
the Local Group ratio. 

5.4. Metals-to-dust ratios 

Metal column densities {Nh,x) from s oft X-ray absorption in 
excess of the Galactic absorption (e.g. iGalama & Wiiersll2001t 



IWatson et al.l2002HStratta et al.l2004tlStarUng et al.l2005ll2007l 



Gendr eet ah 2006; Butler & Kocevskil [20071: IXVatson et al 



2007 



Evans et al 



ISchadv eT aLl l2oTol: ICampana et al.1 l20'Tot 

,20071) have been detected in most GRB afterglows; a rather sur- 
prising fact, since the absorption becomes harder to detect with 
increasing redshift. The effective X-ray absorption is dominated 
by various species, mostly oxygen and siUcon K-shell absorption 
and in cluding helium and iron L-shell dM orrison & McCammonI 
Il983l) . At higher redshifts (i.e. z > 2.5) the soft X-ray absorp- 
tion becomes increasingly domi nated by the heavier elements 
such as iron, silicon and sulphur (Kallman et al.l2003tlGou et al.l 
2005; Butler 2007; Hurkett et al. 2008). It is important to note 
that the X-ray column densities are typically reported in units of 
equivalent hydrogen column density assuming solar abundances. 
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Fig. 6. Metals-to-dust ratio (A^h.xMv) distribution of our spec- 
troscopic GRB sample in Galactic units. The solid line show the 
cases where both A^h.x and Ay are significantly measured. The 
upper and lower limits to the metals-to-dust ratios are shown by 
the dotted histogram and the limits in each bin are marked by 
left or right arrows. 



The upper limits and measured A^h.x values of our spectroscopic 
GRB sample are given in Table |5] Significant excess X-ray ab- 
sorption was detected for 21 GRBs in our sample. The X-ray 
column densities are effectively metal column densities in both 
gas and dust, and would have to be corrected by the metallicity 
to obtain the gas column density. However, this is not enough. 
Due to a lack of information on the precise ionization state and 
structure of the gas, the absorption is invariably fit with a neu- 
tral gas model. We are well aware that the gas surrounding the 
GRB that causes the majority of the absorption is highly ionize d 
dWatson et al.ll2007t ISchadv et al.ll201ll: iProchaska et d1l2007h . 
The metal column densities derived from the fits using neutral 
models are therefore strictly lower limits to the actual metal col- 
umn densities. How close they are to the real value of the to- 
tal column density depends strongly on the ionization state of 
the gas. It is reasonable to assume we are not underestimat- 
ing the column densities by a factor of more than ten, since in 
the cases with the highest column densi ties, Comp ton scatter- 
ing effects would bec ome significant d Watson & Laursea.201 k 
ICampana et alJl20in) . and we would not detect some bursts in 
X-rays. 

In Fig.|5]we plot the dust extinction and soft X-ray-derived 
metal column densities for the sample. It is immediately appar- 
ent that all of the dust extinctions lie far below the expecta- 
tion for the Local Group metals-to-dust ratios. For the Galaxy, 
Nh.x/Av has been mea sured for many sightlines and us ing dif- 
feren t techniques (e.g. Giiver & Qzel 2009; Predehl & SchmittI 
Il995t iGorensteinI Il975h . For the LMC and SMC, the gas- 
to-dust ratio (Nhi/Av) can be converted to the same units 
(Mi.xMv) by correcting for the metallicity. The Nui/Ay relation 
for the MW varies in the range (1.7 - 2 2) x 10^' c i ri'^ ni ag~'. 
Correc ting for metalhcities, the LM C (lFitzpatricklll985b and 
SMC dWeingartner & Draine 20011) value s are consistent with 
this range (see for a review lDrainell2003l) . Therefore, we use a 
canonical value of Nh.x/Av ~ 0.2 x 10^^ cm"^ mag"' for the 
metals-to-dust ratio of the Local Group of galaxies for com- 
parison to our sample. Dai & Kochanek (2009) showed that in 
lensing galaxies, the metals-to-dust ratio is also consistent with 
values in the Local Group. Together, those results hint that there 



may be an almost universal metals-to-dust ratio, at least in galax- 
ies that are somewhat evolved. 

In apparent contrast to this, our sample shows a spread of 
metals-to-dust ratios that covers a range from approximately 
three to thirty times the Local Group value (Fig.|6). The fact that 
none of our objects have a metals-to-dust ratio below the Local 
Group value suggests that the Local Group ratio may at least be a 
minimum. Indeed , our result here, w h ich confirm much previous 
work (e.g ., Gala ma & Wiiersll2001t ISchadv etal.ll2010l 1201 It 
IStratta et al.i2007D , may be consistent with a universal metals-to- 
dust ratio if dust destruction by the GRB and afterglow is respon- 
sible for the lack of dust that leads to these large metals-to-dust 
ratios. Under this hypothesis, the GRB would have to destroy the 
equivalent of between 1 and 10 mag nitudes of extinction (Ay) of 
dust within a f ew pc of the burst ([W atson et al.l 2007h. which 



seems feasible ([Fruchter et~ar ]|200inWaxman & Drainel 



Pernaetal.ll2003h. As p ointe d out bv IFruchter et al.l (120011 



20001 



iWaxman & Draniel(l2000l) . and lPema et al.l (l2003l) . observations 



within the first few seconds of a burst should, for most GRBs, 
show a colour change related to the destruction of the dust. 
Based on these data and the assumption of a universal metal-to- 
dust ratio, such a colour change would be very large; however, 
given that we rely on bright restframe UV emission for detec- 
tion, and it is precisely the bright UV emission that is primarily 
responsible for the dust destruction, it seems unlikely that we 
will be able to detect such a colour change before it is complete. 

It is worth noting here again that neither ionization nor 
metallicity can resolve the dust-poor nature of these absorptions. 
In the former case we know the gas must already be highly ion- 
ized; the dust would not survive in such an environment. As 
to the latter, we are directly observing the metal absorption, 
changes to the overall metallicity therefore have no efi^ect on this 
ratio. 

It is worth noting in this discussion that because our sam- 
ple is almost certainly biased towards low extinction sightlines 
due to the requirement to have an optical spectrum, we may be 
preferentially selecting those sightlines with the highest gas-to- 
dust and metals-to-dust ratios without needing to invoke dust de- 
struction. It is possible that we might simply be missing many 
extinguished GRB afterglows with low metals-to-dust and gas- 
to-dust ratios. However, in this scenario, one would have to ac- 
cept that sightlines in a significant number of GRB hosts ex- 
ist with metals-to-dust ratios tens of times higher than in the 
Local Gro up. Such higher metals-to-dust ratios are a lso found in 
AGN re.g. JMaioUno et alJl200UlWillott et alJl2004h using simi- 
lar techniques to those used here. In those cases, however, most 
of the high X-ray absorbing column is attributed to the AGN and 
is assumed to be essentially dust-free. Sightlines throughout the 
Local Group from the SMC to the Milky Way, on the other hand, 
have metals-to-dust ratios that do not seem to vary by more than 
some tens of percent. 

We examined the possibility of a relationship between the 
metals-to-dust ratio and the line-of-sight extinction (Fig. |7]l. 
While it might be anticipated that GRB hosts that possessed 
high extinction sightlines might be dust-enriched and have low 
metals-to-dust ratios, we see no clear correlation here either. 
Even the highly extinguished events with a 2175 A bump have 
metals-to-dust ratios similar to the rest of the sample. 

We briefly compared the metals-to-dust ratios of our spectro- 
scopic sample with the burst redshift. We see no obvious trend 
of the ratio with redshift except what might be expected from 
the detectability limits of the X-ray absorption increasing with 
redshift and while the limit for Ay essentially decreases with 
redshift. 
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Fig. 7. Restframe Ay and metals-to-dust ratio (A^h.xMv) for our 
GRB sample. The filled circles are for the cases where A^h.x and 
Al/ are > 2cr significant. The open circles represent upper limits 
either obtained for A^h.x or Ay. The filled square corresponds to 
GRB 080607 observed with the Keck telescope. 



5.5. The effect of metallicity 

Metallicities measured in GRBs vary from less than 1/100 
to nearly solar values, and are usually higher than the val- 
ues fou nd in QSO-DLA si ghtlines (see 'Fvnbo et al.' '20064 
[Prochas ka et al.li2007l: iFvnbo et al. 2008c; Ledoux et al. 2009). 
This trend supports the idea that GRBs occur in low mass or 
dwarf galaxies, representing the largest fraction of galaxies in 
the local Universe with large star formation rates. This favours 
the currently accepted scenario of long duration GRB s, the 
collapse of massive, rotating Wolf-Rayet star (Wooslevl ll993t 
iMacFadven & WooslevI Il999l) . The metallicities are obtained 
from Voigt-profile fitting to mildly saturated lines arising from 
non-refractory elements using mid/high resolution spectroscopy. 
In our metallicity sub-sample we obtained metallicities for the 
GRBs from the literature (see Table |6]l. The metallicities based 
on Zn are given for GRB 050401, GRB 070802, GRB 071031 
and GRB 080210. In addition, we derived Si column densi- 
ties, and hence, metallicities, for GRB 060707, GRB 060714, 
GRB 070 1 10, u sing the optically thin limit approx imation (e.g . 
[Petitieaniri998h with oscillator strengths from Mortonl (1 199 ih . 
For these three bursts we used the equivalent width of the weak 
and hence likely only mildly saturated Si ii (1808 A) line to esti- 
mate the Si II column density. For a few more GRBs there is no 
detection of Si ii (1808 A) whereas Si ii (1526 A) is well-detected 
(see F09). We inferred metallicities for eight GRB hosts using 
the restframe equivalent widths of Sin (1526 A). The equiv- 
alent widths for all these GRBs are obtained from F09. The 
metallicities are derived in these ca ses using the metallic ity- 



EWi526 A correlation for QSO-DLAs (iProchaska et al.ll2008l Eq. 
1). Other independent QSO-DLA studies also confirm this tight 
trend between EW,526a and [M/H] (Kaplan et al. 2010); GRB- 
DLAs also exhibit the same trend (IProchaska et al .1 12008b . For 
GRB 080607 the absorption line featur es are highly saturated , 
therefore only lower limits are obtained (IProchaska et al.E009h . 
The metallicity for this burst could be solar or even super solar. 
The reported metallicity here is a lower limit based on Zn. No 
metallicity measurements are available for 60% of our sample, 
mainly due to a lack of high- or medium-resolution spectroscopy 



or a redshift < 2, placing the Lycf line outside the observed wave- 
length range. 

In Fig. [8] we plot metallicities for our GRB sub-sample ver- 
sus the visual dust extinction. Ay. Fig. [8] shows that there may 
be a slight trend between the metallicity and dust extinction. We 
see somewhat higher dust with increasing metallicity, suggest- 
ing that more metal rich systems may be more dusty. However 
we do not claim this cor relation here with any confidence. 

ISchadv etaP dloTol) also found high X-ray-derived metals- 
to-dust ratios and proposed that these ratios could be related to 
the metallicity of the GRB environment. They suggested that 
there might be an anti-correlation between the cold, gas-phase 
metallicity and the metals-to-dust ratios; i.e. metal poor galax- 
ies forming dust less efficiently from their metals than the metal 
rich ones. Fig. [TOl shows the relation using the data in our sam- 
ple. We see no obvious trend between the metals-to-dust ratios 
and the gas-phase metallicities and cannot confirm the relation. 
Fo r the Spitzer Infrared Nearby Galaxies Survey (SINGS) sam- 
ple |Draineeral ] (120071) found that the gai-to-dust ratio was anti- 
correlated with met allicity. Th e anti-correlation was also found 
for dwarf galaxies (iLisenfeld & Ferrara 1998). Such a correla- 
tion would not be terribly surprising since the dust is composed 
of metals. In Fig. |9] we plot the gas-to-dust ratios of our GRB 
sub-sample versus metallicity. There is an anti-correlation be- 
tween the two values; a Spearman rank test for a random origin 
for such an anti-correlation yields a < 1 % probability. So, while 
there is some evidence of an anti-correlation between the gas-to- 
dust ratio and metallicity in GRB host galaxy sightlines, we find 
no evidence for a relationship between the metals-to-dust ratio 
and the metallicity. GRB environments have consistently higher 
metals-to-dust ratios than the Local Group but at the same time 
their metallicities cover the range as observed in the Local Group 
with a large scatter Fig. [TO]in fact shows the large gap between 
the Local Group and GRB hosts in metals-to-dust ratios. 

For most of the GRB afterglows we used Si based metallic- 
ity. In any case, the possible depletion of Si does not affect our 
conclusions above, since if the Si metallicities are affected by 
dust depletion, they will be underestimated. Specifically, high 
metals-to-dust ratios would therefore be higher if corrected for 
dust depletion (Fig.|4|, depletion is likely to be higher in systems 
with lower gas-to-dust ratios, thus making the correlation in Fig. 
|9] potentially stronger, and correcting for possible depletion in 
Fig. [To] is only likely to make the lack of correlation between 
metals-to-dust ratio and metallicity even clearer 

6. Conclusions 

In this work we have presented the results of our analysis of 
41 GRB afterglows from the X-ray to the NIR. This is not only 
the first spectroscopic study of absolute extinction curves in dis- 
tant galaxies, but also the largest extinction sample outside the 
Local Group. The SEDs are almost universally well-described 
by a power-law or broken power-law with absorption in soft X- 
rays and extinction in the optical/UV. Our analysis shows that 
a break occurs between optical and X-ray in about half of all 
SEDs at times typically between a few hours and a day after 
the burst. The mean spectral change excluding one outlier, is 
A/3 = /32 - fix - 0.51, with a very small scatter (0.02), suggesting 
that in future afterglow SED analyses, the spectral change can 
typically be fixed to 0.5. 

63% (26/41) of the bursts in the sample are well-fit with 
an SMC extinction curve, with 27% (11/41) of the afterglows 
found to have no significant extinction. Of the remaining 4, 
3 have a 2175 A bump and one has a very unusual extinction 
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Table 6. GRB sub-sample metallicity measurements taken ei- 
ther from the literature or based on the Si ii lines strength. 
The Sin {A = 1808 A, 1526 A) absorption lines are used to 
obtain metallicity by u s ing th e optically thin approximation 
and the iProchaska et al.l (|2008|) correlation. The metalUcity for 
GRB 080607 is based on Zn element. The metalUcity could be 
solar or even super-solar for this burst (IProchaska et al.lS009h . 



GRB 



[M/H] 



Reference 



050401 

050730 

060115 

060707 

060714 

060906 

060926 

061 HOB 

070110 

070506 

07072 IB 

070802 

071031 

080210 

080605^ 

080607" 

080721 

080905B 



-1.00 ±0.40 
-2.17 ±0.10 
-0.72 ± 0.37 

> -0.40 

> -1.30 
-1.38 ±0.22 
-0.95 ±0.28 
-1.12±0.21 

> -1.40 
-0.31 ±0.14 
-0.67 ± 0.49 
-0.50 ± 0.68 
-1.73 ±0.05 
-0.93 ±0.18 
-0.42 ± 0.29 

> -1.70 

> -1.10 
-1.10±0.26 



Watson et al 




^rochaskaetal^ (2008) correlation 
From Si 11(1808) 
From Si 11(1808) 
Prochaska et al. (2008) correlation 
Prochaska et al. (2008) correlation 
Prochaska et al. (2008) correlation 

From Si 11(1808) 
Prochaska et al. (2008) correlation 
Prochaska et al . (2008) correlation 
Eh'asdottir et al. (2009) 
Ledoux et al. (200^ 
De Ciaetal. (2011) 
Prochaskaetal^ (_2008) correlation 
Prochaska et al. (2009) 
St arling et al. (2009) 
Prochaska et al. (2008) correlation 



" The metallicity for GRB 080605 is based on the Si ll line EW. Being 
at z < 2, H I is not detected for this burst. Therefore it is not included in 
Fig.HandlH 
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Fig. 8. Restframe A V measurements versus GRB metallicity. The 
metallicities obtained from the literature based on Zn are marked 
with red stars. The remaining data points are based on Si el- 
ement metallicities (see Table |6]l. GRB 080607 is shown as an 
open square. The small circles or stars represent metallicity or 
Ay upper limits. 




Fig. 9. Gas-to-dust ratio versus metallicity for GRB sub-sample. 
The metallicities obtained from the literature based on Zn are 
marked with red stars. The remaining data points are based on 
Si (see Table|6]l. GRB 080607 is indicated with a red square. The 
small circles and stars denote metallicity or gas-to-dust ratio up- 
per limits. MW, LMC and SMC environments are denoted with 
blue diamonds from right to left. 
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Fig. 10. Metals-to-dust ratio versus metallicity for the GRB sub- 
sample. The metallicities based on Zn are marked with red stars. 
The black circles represent Si-derived metallicity (see Table |6]l. 
GRB 080607 is marked by a red square. The small circles and 
stars illustrate metallicity or metals-to-dust ratio upper limits. 
MW, LMC and SMC environments are denoted with diamonds 
from light to lef t. The blue triangles coiTespond to four data 
points taken from lSchadv et al . ( 2010). The blue dashed line rep- 
resents the best fit t o the three Local G roup environments and 
the four GRBs from lSchadv et al.1 ( |2010|) . The blue triangle with 
dotted error bars is GRB 050730 which is also included in both 
ISchadv et al.l (1201 Ol) and our sample. 



curve. The bump cases comprise 7% of our spectroscopic sam- 
ple. Most bursts of our sample have moderate or low extinc- 
tion, with Ay < 0.65. However cases with a 2175 A bump all 
have Ay > 1.0. The predominance of the bump in the extinc- 
tion curves of high-extinction sightlines allows us to claim that 
since our sample is heavily biased against high extinction sight- 
lines, missing between a quarter and a half of all bursts, the per- 
centages of extinction curves listed above are not at all repre- 



sentative of GRB hosts as a whole. Previous conclusions sug- 
gesting that the bump is a rare feature in extragalactic sightlines 
are also likely to be incorrect. An important conclusion from 
this work is that the most extinguished systems seem to cor- 
relate with the systems with a 2175 A bump and possibly with 
high metallicity. Both facts point to the most extinguished sys- 
tems being the most evolved, and therefore likely more massive. 
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This is consistent with reports on work by iPerlev et"aD (l2011ah 
based on Spitzer data suggesting that dark bursts are more likely 
to have red, dusty host galaxies. Since a very substantial frac- 
tion (~ 25 - 42%) of GRB afterglows are dark and likely extin- 
guished (F09: iGreiner et al.ll201 lb . and all four known 2175 A 
systems have restframe Ay > 1, corresponding to >2-3 magni- 
tudes of extinction in the observed BVRI bands, it is premature 
to claim unequivocally that GRBs reside in small, sub-luminous, 
blue galaxies, or even that they follow star-formation. 

The gas-to-dust ratios of GRBs are typically higher than the 
SMC. Correcting for gas-phase metallicity does not provide an 
explanation for these high ratios, and we find that gas-to-dust 
ratios from optical/UV measures are likely to be dominated by 
ionization and dust-destruction effects by the GRB, and not in- 
dicative of the intrinsic ratios of the GRB host galaxies. Soft 
X-ray-derived metals-to-dust ratios are 3-30 times the Local 
Group values. This is not due to ionization or metallicity effects. 
GRB hosts may however be consistent with a universal metals- 
to-dust ratio if dust destruction by the GRB is responsible for 
the apparent paucity of dust. We find no correlation between the 
metals-to-dust ratios in GRBs and the cold, gas-phase metallic- 
ity, the redshift or the extinction. 
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Appendix A: Details on the spectroscopic GRB 
afterglow sample 

A.1. GRB 050401 

The optical spectrum of the afterglow o f GRB 050401 was 
secured with the VLT/FORS2 instrument (IWatsonet all 120061: 
F09). The NIR ph otometry in the /, H and K bands from 
I Watson et al.l (l2006h is interpolated to the time of the SED us- 
ing the R-hwA decay rate {a - 0.86 + 0.1). We find that the 
SED of the afterglow can be well reproduced by a broken power- 
law with substantial amount of SMC-type dust extinction with a 
value ofAy- 0.65 + 0.04 mag. The best-fit parameters are given 
in TableS 

GRB 050 401 at z ^ 2.8983 is cate gorized as a "dark burst" 
as defi ned bv lJakobsson et all (l2004ah . Previouslv I Watson et alJ 
(l2006h found the afterglow SED is compatible with a broken 
power-law model and SMC type dust reddening with Ay - 
0.62 + 0.06, in agreement with our results. This burst has re- 
markably one of the hi ghest Hi column densities with log 
A^Hi/cm-^ = 22.60 + 0.3 dWatson et alj|2006h . 



A.2. GRB 050730 

The optical spectrum of the afterglow of GRB 050730 (z - 
3.9693) presented in this paper was carried out with the 
VLT/FORS2 (F09). The afterglow SED is constructed at 4.1 hr 
after the burst. The SED can be explained by a broken power- 
law, consistent with a ~ 0.5, plus a small amount of SMC 
type extinction with Ay = 0.12 + 0.02 mag. 

Optical spectra of GRB 050730 were also acquired by 
the Magellan Inamori Kyocera Echelle (MIKE) echelle spec- 
trograph on the Magellan Telescopes (IChen et al.l l2005h . the 
Intermediate-dispersion Spectroscopic and Im aging System 
(ISIS ) on the William Herschel Telescope (WHT: IStarling et alj 
2005h and VLT UV- visual Echelle Spectrograph (UVES; 



D'Elia et aUlIOOTl: iLedoux et al.ll2009l) . The derived H i column 



density i s log jVHi /cm'^ = 22 1 + .1 dStarhng et al.l 120051: 



iD'Eha e t al. 2007; Ledoux et al]l2009l) . A foreground QSO is 



also detected in the 2D s pectrum at an imp act parameter of 20 
arcsec (F09). Previously (lOates et al"]|2009l) derived Ay - 0.15 
from a fit of UVOT and X-ray data. iKann etaP (1201 Oh find 
Ay = 0.10 + 0.02 from their optical/NIR SED fitting. Both re- 
sults are consistent with our fi ndings. Starli ng et al. (2005|) find 
Ay - 0.01 from spectral fitting. ISchadv et al.l (l2010l) find that no 
extinction model is capable to fit of the NIR through X-ray SED, 
and derive Av/^0.16-0.23 for the SMC model using single and 
broken power-laws. 

A.3. GRB 050824 

The optical spectrum of the afterglow of GRB 050824 
(z = 0.8278) was ta ken with the VLT/FORS2 instrument 
dSoUerman et al.l l2007t F09) GRB 050824 is an X-ray flash 



ISchadv et all (l2007bl) find SMC dust fits the data best, but 
they are not able to discriminate among the single and bro- 
ken po wer-law models. For both power-law models Schady et al] 
(1201 Oh find Ay ^ 0.12-0. 16, in excellent agreement with 
our values. iKann et all (1201 Oh also find that SMC model well 
fits the optical/NIR SED with an insignificant amount of dust 
(Ay = 0.14 + 0.13). 

A.4. GRB 06011 5 

The optical spectrum and the /?-band photometry of the af- 
terglow of GRB ( )60115 (z = 3.5328) we re carried out with 
the VLT/FORSl dPiranomonte et al.ll2006bl: F09). The SED is 
generated at fo + 8.9 hr after the burst. The afterglow SED of 
GRB 0601 15 is described by a single power-law and SMC type 
extinction with a very small amount of dust (Ay = 0.10 + 
0.02 mag). The optical spectrum has very low S/N and this can 
be seen clearly in the extinction curve (see Fig. lA.2) . F09 derived 
an Hi column density of log A^Hi/cm"- = 21.5 +0.1. 

A.5. GRB 06051 2 

The spectrum of the af terglow of GRB 0605 12 was carried out 
with the VLT/FORS2 dStarling et al.ll2006at F09). For this af- 
terglow there is uncertainty about the redshift in the GCN cir- 
culars. The redshift of z = 2.1 comes from a single absorp- 
tion line detected in the Telescopio Nazionale Galileo (TNG; 
IStariing et ani2006bh and the FORSl spectra (F09). This red- 
shift is also consistent with the spectral break detecte d by 
UVOT dde PasQuale & Cummingsll2006HOates et"ani2009h . The 
/?-band after glow light curve is c onstru cted by using magnitudes 
given by .MundeU et al.l d2006h ; ICenko (2006a) ; Milne (.20061) ; 
ICenko & Baumgartner (l2006l) . The A'^-band photometry is in- 
te rpolated at the time of the SED by using observations given 
in iHeartv et al.l d2006l) ; iTanaka et alTd2006l). The 7-band obser - 
vations of the afterglow are taken from Sharapov et al.l (l2006l) . 
We reduced the UVOT data and the afterglow is clearly detected 
in the v, b and u bands. The interpolated lightcurve for each 
band is obtained by using the /?-band lightcurve decay rate of 
Q- = 1.05+0.1. The SED of the afterglow of GRB 060512 fits 
well with a broken power-law and no dust reddening (see Table 
|3]l. We estimated the 2cr upper limit for extinction in the rest- 
frame y-band (< 0.08 mag). 

Previously ISchadv et al.l d2010h find that no dust model is 
capable of fitting the NIR through X-ray data and derive high 
extinction using SMC mode l (Ay ^ 0.47 - 0.66). The extinction 
corrected data presented in ISchadv et al.l (1201 Ol) do not match 
with their given intrinsic slope of the afterglow SED. 



A.6. GRB 06061 4 

The optical spectrum of the af terglow of G RB 060614 was 
obtained with the VLT/FORS2 dDella VaUe et al.. ,2006; F09). 
GRB 060614 has a redshift of z = 0.1257, flie lowest 
redshift in our flux calibrated GRB sample. The redshift 
is estim ated from the emissio n lines from the GRB host 
galaxy (iDella VaUe et alJ l2006h . GRB 060614 is a nearby, 
long duration GRB but not accompanied by a bright super- 



(XRF) dSoUerman et al.ir2007h . The optical and X-ray spectrum nova (iFvnbo et al.l20064lDeUa VaUe et al.l2006l:lGal-Yam et all 



are normalized to the level of acquisition image photometry 
taken at 9.5 hr after the burst trigger The SED of the afterglow of 
GRB 050824 is well explained by a broken power-law and SMC 
extinction with Ay = 0.14 + 0.04 mag. 



120061: iGeh rels et al.'2006h. The UVOT data w ere previously pub- 
lished in.Xuetal. (2009!) and iMangaiioet al. (2007). We re- 
reduced the UVOT data for v, b, u, uvwl, uvml and calibrated the 
afterglow fluxes in these bands using the photometric technique 
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given in lPoole et all (l2008h . The optical to the X-ray SED of the 
afterglow of GRB 060614 is well modeled with a broken power- 
law and low SMC type extinction with A ^ = 0.1+ 0.03 mag. 

iMangano et alJ(l2007h find that SMC dust well fits the optical 
through X-ray SEP , with Ay = 0.05 + 0.02, values comparable 
to ours. lKann et al.l (1201 lb find SMC dust model fits the optical 
data, and even higher extinct ion is obtained (Ay = 0.28 + 0.07). 
The SED is not provided in iKann et aP (1201 ih . therefore a de- 
tailed comparison could not be made. 

A.7. GRB 060707 

The optical spectrum of the afterglow of GRB 060707 (z - 
3.4240) toget her with j?-band photo metry are taken with the 
VLT/FORSl (iJakobsson et alJl2006al: F09). The afterglow SED 
is constructed at 34.4 hr after the burst trigger The optical to 
X-ray SED of the afterglow of GRB 060707 can be nicely re- 
produced with a single power-law and SMC type dust extinction 
(Ay - 0.08 + 0.02 mag ). The derived H i colum n density is log 
Nui/cm-^ = 21.0 + 0.2 (IJakobsson et al.ll2006al) . 

A.8. GRB 060708 

The optical spectrum and /?-band photometry of the after- 
glow of GRB (360708 are obtained with the VLT/FORS2 
(IJakobsson et al.ll2006cl : F09). A precise redshift of the burst is 
not obtained due to the low S/N of the optical spectrum (F09). 
The reported redshift of z = 1 .92 wa s based on sp ectral break 
seen in UVOT data due to Lyor break (ISchadv & Mo retti 2006). 
Together with the X-ray and optical spectrum we reduced the 
UVOT data in all bands obtained from the Swift science archive 
facility. The data were reduced using the st andard UVOT pho - 
tometric calibration technique described in iPoole et all (|2008|) . 
The afterglow is clearly detected in the UVOT v, b and u photo- 
metric bands. The optical to X-ray SED of the afterglow is well 
fit with a broken power-law together with SMC type dust extinc- 
tion. Ay = 0.13 ±0.02 mag. 

A.9. GRB 06071 4 

The optical spectrum of the afterglow of GRB 060714 [z - 
2.7108) was carried out with the VLT/FORSl. The TJ-band 
lightcurve is obtained from the circulars spanning from 0.04 to 
4 days after the burst. The / and 7-band photometric magnitudes 
were obtained from the CTIOA Novel Dual Imaging CAMera 
(ANDICAM) instrument (ICobbll2006l) . The UVOT data of the 
afterglow is reduced and the source is well detected in the v and b 
bands. We modeled the intrinsic SED of the afterglow and found 
that the afterglow can be well explained by a broken power-law 
together with SMC type extinction with restframe extinction of 
Ay = 0.22 ±0.02 mag. 

Previouslv lJakobsson et al.l (l2006al) obtained the redshift and 
the Hi column density of log A^Hi/cm"^ = 21.8 ± 0.1 from 
the optical afterglow spectrum. Moreover extended Lya emis- 
sion in the center of the DLA trough is clearly detected (F09). 
ISchadv et al.l (l2010l) find that no dust model is able to fit the 
optical/X-ray data. The least bad is the SMC-type dust and 
power-law with dust extinction Ay = 0.46 ± 0.17, compa ra- 
ble to our results. The afterglow SED in lSchadv et al.l (l2010l) is 
constructed from the Rvb band and X-ray data. While b band is 
clearly affected by the Lya absorption at z = 2.7108 and might 
causes high amount of extinction. 



A.10. GRB 060729 

The optical spectrum of the afterglow of GRB 060729 at z = 
0.5428 is obtained with the VLT/FORS2 ( iThone et alJl2006H: 
F09). The acquisition image photometry was carried out using 
the /-band filter as reported in Table |T] The UVO T data of the 
afterglow of GRB 060729 is previously published in lGrupe et al.l 
(2QQ7). The Swift UVOT data of the afterglow of GRB 060729 
have been re-reduced and the source is clearly detected in the 
V, b, u, uvwl and uvm2 filters. The optical spectrum is not of 
good quality and has broad undulations, therefore, the true shape 
of the underlying continuum cannot be seen from the spectrum. 
The optical spectrum has very low S/N and is affected by the 
high airmass at the time of the observation (F09). The afterglow 
SED can be fitted by a single power-law and SMC type dust 
extinction, and resulted in a very small amount of dust (Ay = 
0.07 ±0.02 mag). 

ISchadv et aP (1201 Ol) find that SMC/power-law is the best 
model, resulting in low amount of extinction with Ay = 0.03 ± 
0.01 and consistent with our value. 



A.11. GRB060904B 

The optical spectrum of the afterglow of GRB 060904B 
z = 0.7029) was o btained with the VLT/FORSl instrument 
Fugazza et al.l l2006t F09). In our sample GRB 060904B has 
the highest amount of foreground extinction and absorption (see 
Table [TJ. The SED of the afterglow is constructed at 5.1 hr af- 
ter the burst. The SED fits well with a single power-law and 
SMC type extinction with extinction of Ay = 0.31 ± 0.02 mag. 
The SED seems like a broken power-law case with softer X-ray 
slope. However, due to large uncertainties on the X-ray spectrum 
the fitting routine cannot minimize both models accurately and 
prefers a single power-law. 

Previouslv lSchadv et al.l (1201 Ol) find no good fit for any dust 
model but find that single power-l aw is less likely, su ggesting 
low extinction with Ay = 0.12+° °^. lKann et alj ( 1201 Ol) find that 
their optical-NlR SED is consistent with no d ust reddening usin g 
any dust model. The SED is not provided in iKann et alJ ( 1201 Ot) . 
therefore a detailed comparison could not be made. 

A. 12. GRB 060906 

The optical spectrum of the afterglow of GRB 060906 [z = 
3.6856) was obtained with the VLT/FORSl spectrograph 
( IJakobsson et al.l2006al: F09). The /-band observation was taken 
at the New Mexico Skies Observatory (Torii 2006) and corrected 
for the time using the /?-band lightcurve power-law w ith shallow 
decay index a = 0.56 ± 0.07 dCenko et al.ll2006l:ICenko,.2006bl) . 
The z' band observ ation is taken from the lightcurve presented in 
ICenko et al.1 (|2009|) . The optical to X-ray SED can be fitted with 
a broken power-law with no dust extinction. The derived 2cr ex- 
tinction upper limit is < 0.09 mag. The X-ray to optical slope 
differ by A/3 ~ 0.5 which corres ponds to the ch ange in slope 
caused by the cooling frequency (ISari et aLlll998l) lying within 
the observed frequency range. 

The derived H i column density is l og NuJcm^^ - 21.8 5 ± 
0.1 (IJakobsson et al.ll2006al) . Previouslv ICenko et al.l (l2009l) fit- 
ted the SE P and found dust extinction with Ay = 0.20+° j2. 
iKann et aP ( 12010) fit the optical afterglow SED and find in- 
significant dust extinction Ay = 0.05 ± 0.05, in agreement with 
our results. 
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A. 13. GRB 060926 

The optical afterglow of GR B 060926 (z = 3.2086) w as ob- 
served with the VLT/FORSl dPiranomonte et al.ll2006i F09). 
The Lya emission from the host galaxy of the GRB is detected 
in the trough of the DLA (F09). T he estimated H i column den- 
sity is log A^Hi/cm"^ = 22.7 + 0.1 (iJakobsson et al.ll2006ah . The 
optical to X-ray SED is normalized to the acquisition image pho- 
tometry at 7.7 hr after the burst. The optical to X-ray SED is well 
described by a single power-law and SMC type dust with a mod- 
erate extinction of Ay = 0.31 + 0.02 mag. 



A.16. GRB 061 021 

The spectrum of the G RB 061021 was se cured with the 
VLT/FORSl spectrograph (iThone et al.ll2006al: F09). The red- 
shift of the burst (i.e. z - 0.3463) is obtained from Mg ii ab- 
sorption lines (Hjorth et al, in prep). This is the second low- 
est redshift of our flux calibrated GRB afterglow sample. The 
UVOT photometric data have been reduced and the afterglow is 
well detected in the v, b, u, uvwl and uvm2 filters (see Fig. lA.lb . 
The intrinsic SED of the afterglow can be explained by a bro- 
ken powe r -law c onsistent with a cooling break as explained by 
ISari et aL l (fT998h and with no dust extinction. The estimated 2cr 
extinction upper limit is < 0.1 mag. 



A.I 4. GRB 060927 

The optical spectrum of the afte rglow of GRB is car ried out with 
the VLT/FORSl spectrograph (iRuiz-Velasco et al.ir2007 ; F09). 
GRB 060927 with a redshift of z = 5.4636 is the second highest 
redshift GRB in our spectroscopic G RB sample. The redshift i s 
based on a single Si ii absorption line (iRuiz-Velasco et al l2007h . 
The K, J and / band observations of the aftergl ow were carried 
out with VLT filters (iRuiz-Velasco et al.ll2007h . We took NIR 
photometry and spectroscopy to construct the composite SED at 
12.5 hr after the burst trigger. The SED of the afterglow can be 
fitted with a single power-law without any dust extinction. The 
estimated 2cr upper limit for the Ay is < 0.12 mag. 

The optical spectrum of t he afterglow is previously published 
in lRuiz-Velasco et"an (l2007h rep orting the H i column density of 
log A^Hi/cm-2 = 22.50 + 0.15. iKann et alj ( I20T0I) fit the NIR 
SED of the afterglow and found < 2.5cr significant dust with 
Ay = 0.21 + 0.08 for SMC dust. Due to the large uncertainties 
on the NIR data, no model could be able to fit the NIR segment 
alone. 



A.17. GRB 061 11 OA 

The optical spectrum of the afterglow of G RB 061 1 lOA at z = 
0.7578 was obtained with the VLT/FORSl jFvnbo et al. 2003; 
F09). The intrinsic SED of the afterglow is scaled to the acquisi- 
tion image photometry. The afterglow spectrum of this burst has 
low S/N. The optical to X-ray SED can be modeled with a single 
power-law with no dust reddening. We derived the 2cr extinction 
upper limit of Ay < 0.1 mag. 

A.I 8. GRB 061 11 OB 

The VLT/FORSl caiTied out the optical o bservation of the af- 
terglow of GRB 061 HOB at z = 3.4344 ( iFvnbo et al.ll2006al: 
F09). The afterglow SED is constructed at to + 2.5 hr epoch 
scaling with the /?-band flux. The estimated Hi column den- 
sity is log Nui/cm-^ = 22.35 + 0.1 (F09). The intrinsic SED 
of the afterglow can be explained by a broken power-law all the 
way from optical to the X-ray with SMC type extinction with 
Ay = 0.22 ±0.03 mag. 



A.I 9. GRB 0701 10 



A.I 5. GRB 061007 

GRB 06 1007 is an extremely bright burst detected by Swift ac- 
companied by a very luminous afterglow and simila r decay rate 
in the X-ray and optical bands ( ISchadv et al.ll2007al) . The spec- 
trum of the optical afterglow of GRB 061 007 at z = 1.2 622 was 
obtained with the VLT/FORSl CJakobsson et alj|2006bl: Paper). 
Th e /-band photometry is obtained from the lightcurve presented 
by iMundell et all (l2007h . We reduced the UVOT data of the 
GRB afterglow and detected source in the v, b and u photo- 
metric bands. The SED of the afterglow is nicely fitted with a 
single power-law together with SMC extinction with a value of 
Av = 0.35 +0.03 mag. 



lOates et alj ( I2009I) and lSchadv elal] ( 1201 Ol) fin d high values 
of exti nction varying from Ay ~ 0.66 - 0.75, while Schadv et al 



201c ) prefer LMC dust model. The results of Mundell et al 
(I2OO7I) with Ay = 0.48 + 0.19 and iKann et all (l2010l) wkh 
Ay - 0.48 ± 0.10 are comparable with our values. There is an 
indication of a 2175 A bump in the blue end of the optical spec- 
trum. The central wavelength of the bump (i.e. /lobs ~ 4520 A) 
does not match with the redshift of the GRB. An intervening ab- 
sorber at redshift z = 1.066 is also seen in the spectrum with 
strong Mg II and Can lines. This absorption feature could be 
associated with the intervening absorber. Apart from the optical 
spectrum, no other sign of the intervening absorber is seen at this 
redshift. 



The optical spectrum of GRB 0701 1 (z - 2.3521) was obtained 
with the VLT/F0RS2 spectrograph (iJaunsen et al.ll20b7b- F09). 
The UVOT data of this burst have been reduced b y using UVOT 
calibration technique given bv lPoole et al.l(l2008h . The afterglow 
is detected only in the v and b bands. The optical to X-ray SED 
of the afterglow is fitted with a broken pow er-law indicating a 
cooling break as sugg ested bv lSariet al.l ([T998l. We find no dust 
extinction in this burst. The derived 2cr extinction upper limit is 
Ay < 0.11 mag. 

The estimated Hi c olumn density i s log A ^h i/cm~^ = 21.7 + 
0.1 (F09). Previously ISchadv et all (1201 Ol) find that broken 
power-law well fits the data, but they could not distinguish be- 
tween the MW, LMC and SMC dust models, finding dust extinc- 



tion of Ay = 0.23+ 



A.20. GRB 0701 25 



with SMC model. 



The optical spectrum of the afterglow of GRB 070125 (z = 
1.5471) was carried out with the VLT/FORS2 and previously 
published in F09. The broad undulations in the optical spec- 
trum are due to systematics in the flux calibration. The red- 
shift of the GRB is determined from the data o btained with 
Gemi ni North Telescope equipped with GMOS (ICenko et al.l 
I2008h . The Rc and / band observations have been collected 
by the MITSuME Telescope ( Yoshida et al. 2007). The /, H 
and Ks band observations of the afterglow have been obtained 
with the robotic Peters Automatic Infrared Imaging Telescope 
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(PAIRITEL; iBloom et alJ l2007bh and scaled to Rc-hand flux 
level using the decay rate from a decent /?-band lightcurve 
obtained from the circulars. The /?-band lightcurve exhibits a 
plateau phase and a break to a steeper decay with decay in- 
dices ffi - 1.4 + 0.1 and q'2 - 2.55 + 0.07 and break time at 
~ 4 days after the burst. The optical spectrum is normalized to 
7?c-band flux level. The X-ray to optical/NIR SED can be fitted 
with a broken power-law and SMC-type linear extinction with 

Av = 0.27 + 0. 03 mag. 

Previously iKann et al ] 1I2OIOI) carried out a fit to the opti- 
cal/NIR data and found fliat SMC dust best fits the data with 
Uttle amount of dust (Ay = 0.1 1 + 0.04). We find that the errors 
on the photometry are large and the NIR/optical X-ray joint fit 
suggests a cooling break in the SED and moderate extinction. 



SED of the afterglow is modeled and we find a broken power- 
law to be a reasonable fit to the data. In addition to the broken 
power-law the data require dust using SMC type dust grains with 
= 0.44 + 0.05 mag. 

A.24. GRB 070611 

The optical spectrum of the afterglow of GRB 0706 11 {z = 
2.0394) was obtained with the VLT/FORS2 (iThone et al.l2007al: 
F09). The derived neutral hydrogen column density of the DLA 
is log Nm/cm-^ = 21.3 + 0.2 (F09). The SED is scaled to the 
time of the /?-band photometry at 7.7 hr The X-ray to optical 
SED of the afterglow is defined by a single power-law and SMC 
type dust extinction {Ay - 0.06 ± 0.02 mag). 



A.21. GRB 0701 29 



A.25. GRB 070721 B 



The optical spectrum of the afterglow of GRB 070129 was ob- 
tained with the VLT/FORS2 spectrograph together with the R- 
band imaging (F09). We collected the /-band detection of the 
aftergl ow from Michigan-D artmouth-MIT (MDM) Hiltner tele- 
scope dHalpem et al.l 120071) . which is close to the time of the 
SED. The redshift (z = 2.338) reported here was estimated from 
the host galaxy [O ill] emission line (Bo Milvang-Jensen private 
comm.). The intrinsic SED of the afterglow is scaled to the level 
of the /?-band observation. The X-ray to optical SED of the af- 
terglow is reproduced well with a broken power-law with SMC 
type extinction with a typical value of Ay - 0.30 ± 0.02 mag. 

A.22. GRB 07031 8 

The optical spectrum of the afterglow of GRB 0703 18 (z - 
0.8397) was carried out with the VLT/FORS 1 spectrograph 
dJaunsen et al.l 12007 at F09). The spectrum shows an unusual 
sharp break in the optical spectrum around A^bs ~ 5000 A. A 
spectrum simultaneously obtaine d at the M agellan telescope, 
showed the same optical break ( IChen et al.l l2007i) . The Swift 
UVOT data have been reduced together with the NIR data from 
ANDICAM. The UVOT slope is consistent with the downturn in 
the optical spectrum. The SED of the afterglow is generated at 
?o + 1 6 .7 hr after the burst. There is clear evidence for dust attenu- 
ation in the NIR to X-ray afterglow SED. Due to the sharp break 
in the optical spectrum and peculiar dust reddening the data are 
not consistent with SMC-origin extinction curve. We fit the SED 
with FM-dust induced model and found that the SED fits with 
a sharp break in the optical spectrum. The extinction curve is 
very unusual due to the sharp optical break. The GRB is not in- 
cluded in our extinction sample due to non-consistency with the 
synchrotron model. This sharp break could be due to destruction 
of dust grains by the GRB. A more detailed and comprehensive 
analysis of the SED will be g i ven in Watson et al. (in prep). 

Previously ISchadv et all ( 1201 Ol) implemented UVOT-XRT 
joint fit and found that a broken power-law is a slightly better 
fit to the data with dust extinction of Ay = 0.59+° °^. 



The optical spectrum of the afterglow of GRB070721B (z = 
3.6298) is carried o ut with the VLT/FORS2 spectrograph 
(iMalesani et al.l 120071 : F09). The estimated Hi column density 
is log Nm/cmr^ = 21.5 ± 0.2 (F09). The SED is scaled to the 
flux level of the /?-band photometry. The intrinsic SED of the 
afterglow is nicely reproduced by using a broken power-law and 
SMC-type extinction with Ay = 0.2 ± 0.02 mag (see Fig. lA.lb . 

A.26. GRB 070802 

The optical observations of GRB 070802 (z = 2.4541) were ob- 
tained with the VLT/FORS2 spectrograph. The 2175 A dust ex- 
tinction featur e is clearly seen in the optical spectrum of the af- 
terglow (F09; ' Elfasdottir et al. 20091). The bu r st is ca tegorized 
as dark burst as defined by Jakobsson et al.l (l2004al) . The J, 
H and K band observa tions are taken at tn + 2.0 hr from the 
lightcurves presented in lKriihler et al.l (|2008|) for each band. We 
reduced UVOT da ta using the calibration techniques defined by 
iPoole et aP (l2008l) . The afterglow is well detected in the v and 
b bands (see Fig. lA.lb . Due to the presence of a strong 2175 A 
bump feature, the SMC extinction curve gives the worst fit to 
the data. The observed bump is shallower than the MW, there- 
fore, CCM is also not a good fit to the data. We fitted the SED 
using a FM extinction model with both single and broken power- 
laws. The intrinsic SED is well described with a single power- 
law and a large amount of dust observed in rest-frame V-band 
(Av = 1.19 ± 0.15 mag) with total-to-selective extinction of 
Ry - 2.81 ±0.68. The complete log of the best fit FM parameters 
is given in Table |4] 

The neutral hydrogen column density of log 
A^Hi/cm"^ - 21.5 ± 0.2 is d etermined from the optical 
spectrum (lElfasdottir et al.l l2009l) . The rest-frame extinction 
curve of the afterglow is shown in Fig. IA.2I T he NIR t o X-ra y 



SED has been pre v iously fitted by E lfasd ottir et al.l ( 2009 1: 
iKruhler et all (l2008l) : ICenko et al.l (l2009h : iGreiner etTll (I2OI iL 
finding high dust extinction (Ay > 1 .0) with a single power- law, 
similar to our results. 



A.23. GRB 070506 

The optical spectrum of the afterglow of GRB 070506 
(z = 2.3090) was taken with the VLT/FORS2 instrument 
( Thone et alj|2007bl: F09). The estimated Hi column density is 
log A^Hi/cm"^ - 22.0 + 0.3 (F09). The optical spectrum red- 
wards of /lobs ~ 7000 A is affected by strong fringing, there- 
fore, this part is not included in our SED analysis. The intrinsic 



A.27. GRB 071 020 

There was no PC mode data available near the time of the optical 
spectrum of GRB 071020 (z = 2.1462), therefore, we reduced 
early and late time PC mode X-ray data for this burst using the 
procedure described in 32.41 We further checked that the X-ray 
spectral slope for early and late time are similar. The optical 
spectrum of the afterglow was obtained with the VLT/F0RS2 
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instrument (Ijakobsson et al.ll2007d: F09). The derived Hi col- 
umn density is log A^Hi/cm"^ < 20.30 (Nardini et al. in prep; 
F09). In our spectroscopic sample, this is the lowest column 
density obtained from the optical spectrum. The Lya is in the 
very blue end of the spectrum and is not seen clearly (see F09). 
The/?, J and K band observations are taken from Nardini et al. 
(in prep). The //-ban d photometry is carried out with PAIRITEL 
(Bloom et al.il2007al) and scaled to the time of the SED. A com- 
prehensive analysis of the optical spectrum of GRB 071020 will 
be discussed in Nardini et al. (in prep). The SED of the afterglow 
of GRB 071020 can be modeled using a broken power-law with 
SMC type extinction with Ay = 0.40 + 0.04 mag. The change in 
slope is consistent with A/3 = . 5 (see Table[3l). 

The best fit of iKann et al] ilOldl) implies SMC-type dust 
with Av - 0.28 + 0.09, comparable to our results within un- 
certainties. 



A.28. GRB 071 031 

The acquisition image photometry and the spectroscopy of the 
afterglow of GRB 071031 ( z = 2.6918) were c arried out with 
the VLT/F0RS2 instrument (iLedoux et al.l2007t F09). The z', 
J, H and Kg band ph otometry have been obtained from GROND 
jKriihler et al.l2009h . We also compared acquisitio n camera pho- 
tometr y to the r' band photometry presented in iKriihler et al.l 
(l2009l) " and taken at similar time, finding consistent values. The 
obtained colu mn density of the G RB-DLA is log A^ei/cm"^ = 
22.15 + 0.05 (iLedoux et al.ll2009h . The SED of the afterglow of 
GRB 071031 fits well with a broken power-law and no dust red- 
dening. The estimated 2cr reddening upper limit in the V-band is 
< 0.07 mag. 

Previouslv iKann et alj ( 1201 (]|) modeled the optical-NIR data 
and found that SMC model well fits the data with insignificant 
amount of dust Ay - 0.14 + 0.13. The SED pres ented in this 
paper clearly indicate no dust. iGreiner et al. find no dust 

extinction from the GROND-XRT joint fit, consistent with our 
results. 



A.29. GRB071112C 

The optical spectrum of the afterglow of GRB 071 1 12C at z - 
0.8227 was carried out with the VLT/FO RS2 instrument together 
with the acquisition image photometry (Ijakobsson et al.ll2007bt 
F09). The only av ailable H band photometry has been taken with 
the MAGNUM (Minezaki et al. 2007). The /?-band Ughtcurve 
is retrieved from the GRBlog and fitted with a decay index of 
a - 0.92 + 0.06. We used the /?-band lightcurve decay rate to 
obtain the NIR photometry at 10 hr after the burst. The SED and 
optical/NlR observations o f the a fterglow have been previously 
discussed in lUehara et"ai] (l20Toh . The intrinsic SED of the af- 
terglow is modeled well with a broken power-law and presents 
no dust extinction. The derived upper limit for dust reddening is 

< 0.08 mag. 

Previouslv iKann et al.l (|2010') studied the optical-NIR after- 
glow SED and found that SMC dust best fit the data with an 
insignificant amount of dust (Ay = 0.23 + 0.21). On the con- 
trary, NIR to X-ray afterglow SED studied in this paper suggests 
no dust extinction (see Fig. lA.ll) . 



l2007at F09). The estimated redshift of the afterglow comes from 
[On] emission line. The optical and X-ray spectrum are nor- 
malized to the R-hmd acquisition camera photometry at 9 hr 
after the burst trigger The NIR to X-ray data fit well with a 
single power-law and SMC type with dust reddening of Ay = 
0.26 + 0.02 mag. 



A.31. GRB 080210 

The optical spectrum of the afterglow of GRB 080210 at 
z = 2.6419 w as carried out using the VLT/F0RS2 instrument 
(Ijakobsson et al . 2008b; De Cia et al. 201 1; F09). The early and 
late time /?-band ph otometry has been acquired with FORS2 (see 
iDe Cia et aTll201 ll) . We selected the /?-band photometry at 1.69 
hr which is near the mid time of the optical spectrum and scaled 
the afterglow SED to that level. The optical spectrum is cor- 
rected f or slit losses and ki ndly provided by Annalisa De Cia 
(see also lDe Cia et alj|201 lb . In addition to the optical spectrum 
the optica l/NIR data have bee n obtained by using the GROND 
telescope ( Greiner et al.ll20lTl) . We obtained the g', r', /', z', J, 
H and band photometry at the time of the SED (Kriihler, pri- 
vate comm.). The and z' data show a spectral plateau and are 
clearly down compared to the rest of the photometry (see Fig. 
lA.ll i. therefore, the GROND photometry is not included in our 
SED fitting. The optical-XRT SED is modeled nicely with a bro- 
ken power-law and SMC-type dust reddening. The dust extinc- 
tion is found to be Ay - 0.33 + 0.03. The change in the spectral 
slope for this afterglow is AyS = 1.14 + 0.13 (see Fig.[TJ. 

The estimated Hi column density for the DLA is of log 



ty 

Nmlcmr^ = 21.90±0.10 (lDe Cia et al.boTlh . lKann et alJ ( l201(]h 
fitted the GROND SED and assumed the /' and z' data as an in- 
dication of 2175 A bump, resulting in some deviation from the 
dust model and give quite large extinctio n {Ay - 0.7). An in- 
dependent analysis of iDe Cia et a also confirms that the 
spectral change is around 1.0 and is inconsistent with the fire- 
ball model. This is the only outlier with twice the spectral slope 
change as compared to the other break frequency cases. The X- 
ray slope for this case is much softer than the other cases (see 
Fig.[T]i. We checked that the X-ray spectrum becomes gradually 
softer (r - (3+1 > 2.0) in the time interval ranging from 1000 to 
10000 s, suggesting that energy injection from the central engine 
is li kely to be de c reasin g at these times. For SM C type extinc- 



tion 



De Cia et al.l (1201 lb found Ay = 0.18 + 0.03. [Greiner et al.l 



A.30. GRB 071 117 

The optical spectrum of the afterglow of G RB 071 117 (z - 
1.3308) was obtained with the VLT/FORSl dJakobsson et alJ 



(1201 ll) emploved a joint GROND-XRT fit to the data with fixed 
Ay6 = 0.5 and find Ay = 0.24 + 0.03. 

A.32. GRB 08031 9B 

The optical spectrum of the afterglow of GRB 0803 19B (z = 
0.9382) was obtained with the VLT/FORS2 spectrograph al- 
most one day after the burst trigger (F09). The GRB is com- 
monly referred to as th e "naked eye burst" due to its peak visual 
magnitude o f 5.3 m ag (iRacusin et al.ll2008l: I Wozniak et aDl2009l; 
'Bloom et al."2009'). The g', r', i' and z' band photometry has 
been obtained from the lightcurves for these filters presented in 
iTanvir etld] (l2010l) . The UVOT data have been reduced and the 
magnitudes obtained for the v, b, u and uvwl filters, where the 
afterglow is clearly detected. We generated the SED at 26 hr after 
the burst and modeled it using power-laws and dust extinction. 
We found the observed data fit very well with a single power-law 
all the way from optical to the X-ray and no dust reddening. The 
upper limit on the amount of dust is estimated to be < 0.1 1 mag. 
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iKann et all ( 1201 Ol) find no dust extinction from optical-NIR 
SED fitting, in excellent agreement with our result. 

A.33. GRB 080520 

The optical spectrum of the afterglow of GRB 080520 (z - 
1.5457) was obtained with the VLT/FORS2 spectrograph at the 
time when the afterglow was very faint. The UVOT data have 
been reduced and the afterglow is only detected in the v band. 
The SED of the afterglow is constructed and the data are well fit- 
ted with a single power-law and SMC-type dust extinction with 
Ay ^ 0.23 +0.02 mag. 

iGreiner et all (1201 ll) fit the XRT through optical/NIR data 
and find weakly constrained high amoun t of extinction (Ay - 
0.53^Q 42). Our results are consistent to lOreiner et aP (1201 ll) 
within Icr. 

A.34. GRB 080605 

The optical spectrum of the afterglow of GRB 080605 (z - 
1 .6403) was carried out with the VLT/FORS2 in different grisms 
dJakobsson et al.l l2008cl: F09). The 2175 A absorption dip is 
clearly seen in two FORS2 grism spectra. Being at redshift z < 2 
the Lya is not detected for this burst due to the sensitivity range 
of the spectrograph. Therefore metal abundances cannot be ob- 
tained for this interesting case. We reduced the UVOT data of 
the afterglow and found that the source is contaminated by a 
nearby object. After subtraction the afterglow is not visible in 
any UVOT bands. The 2175 A absorption dip is very weak and 
narrow. The SED is modeled with the FM extinction law. The 
SED is well described by a single power-law and large amount 
of dust extinction in restframe V band ofAy- 1 -2^0 'j'o mag. The 
extinction curve of the afterglow suggests Ry consistent with 
its value in the MW (see Fig. IA.2b . The metallicity obtained 
from the equivalent width of Si ii (1526 A) suggests a fairly high 
me tallicity for this burst (see Table |6). 

iGreineret aD(l201ll) find that the MW dust model best-fit the 
data with moderate amount of extinction (Ay = 0.47+0.03). The 
SED presented in this paper is not complete and lacking the NIR 
data, which could be in the restframe V-band. This particular 
case with other 2175 A bump cases will be discussed in Zafar et 
al. (in prep). 

A.35. GRB 080607 

The optical spectrum of the afterglow of GRB 080607 [z - 
3.0368) was obtained with th e Keck telescop e and previo usly 
published in Proch a ska et al.l (|2"009); Perle v et"all (1201 Ibl) and 
F09. Because of being not observed with the VLT, this burst is 
not a part of our spectroscopic GRB sample. The derived neu- 
tral hydrogen column density of GRB -DLA is log A^Hi/cm"^ - 
22.70 ± 0.15 (iProchaska et al.1 120091 F09). Due to the occur- 
rence of GRBs in star-forming molecular clouds, molecular 



non-detections (jVreeswiik et al. 


12004; 


Tumhnson et al.l 120071 


ISheffer et al.ll2009t). Fvnbo et al. 


(12006b 


1 tentatively interpreted 



an absorption feature as H2 in their spectrum of GRB 060206. 
For the first time, strong unambiguous H2 and CO molecular ab- 
sorption lines were detected in the optical afterglow spect rum of 
GRB 080607 ( IProchaska et al.llIOM: jPerlev et al.ll201 Ibl) . Most 
of the 2175 A extinction bump is clearly detected in the opti- 
cal spectrum. The remaining part of the bump is covered by 
using the z and / filters. The afterglow is detected in the NIR 



J, H and K, bands with PAIRITEL. The / and V band ob- 
servations have been performed with the Katzman Automatic 
Imaging Telescope (KAIT) and the R, i and z band detection 
are from th e robotic Palomar 60 inch telescope (see Table: 1 
iPerlev et"d![ 201 Ih). This burst cl assifies as dark by the def- 
inition of Jakobsson et al. ] (l2004al) . The GRB afterglow is re- 
markably luminous in optical and NIR wavelengths. Since the 
optical data have strong H2 molecular absorption lines which 
makes it hard to see the underlying continuum, therefore, we 
fit with the bi nned spectrum whe re lines are already taken out 
(see Table: 4: lPerlevet ai]|2011bh . We fit the data with the FM 
parameterization and the best fit parameters are given in Table 
|4| The intrinsic SED of the afterglow is nicely reproduced with 
Ay = 2.31+Jj-^^mag and a very high Ry = 3.75 ± 1.04. The 

2175 A is the widest compared to all other bump cases (see Table 
HJ. The host galaxy of the G RB afterglow is al so extremely red 
in the optical a nd NIR bands dChen et al.ll2010l) . 

Previously Prochaska et al. I (I2OO9I) found Ay X! 3.2 and 
iPerlev et all (l2011bl) estimated A u = 3.3 + 0.4 the optical/NIR 
data fit. Both authors found Ry - 4 which is consistent with our 
result of a shallower extinction curve. 



A.36. GRB 080707 

The optical spectrum of the afterglow of GRB 080707 (z = 
1.2322) is ob tained with the VLT/FORS 1 under bad conditions 
dFvnboet al.ll2 008b: F09). The SED is normalized to the /J-band 
observation taken at 1.1 hr after the burst. The intrinsic SED of 
the afterglow is fitted with a broken power-law and no dust red- 
de ning. The uppe r Umit on the dust content is Ay < 0.13 mag. 

IGreineret all (120 lib implemented the GROND-XRT joint fit 
to the data and find that the data are consistent with no dust with 

f0.14 
-0.08 



an insignificant value of A v = 0.11^° '*' 



A.37. GRB 080721 

The optical spectrum of the afterglow of GRB 080721 at 
z = 2.5914 was ob tained with the VLT/FORS 1 spectrograph 
(IStarling et al.ll2009l: F09). We extracted the UVOT data and the 
afterglow is detected in the v and b band. The /-band photom- 
etr y is obtained at 10.2 hr using the /-band observations given 
bv IStarling et al.l ( l2009l) . The neutral hydrogen column density 
of the DLA is log Nm/cm-^ = 21.6 + 0.1 (F09). The intrin- 
sic SED is scaled to the /?-band acquisition camera photometry. 
The data is well described by using a single power-law and no 
dust reddening. The estimated upper limit for dust reddening is 
<0.12mag. 

Previouslv lStarling et al.l (l2009l) found no cooling break be- 
tween the X-ray and op tical wavelengths. Despite the large er- 
rors on the photometry, iKann et all (1201 Ob fit the optical SED 
and find Ay =0.35 + 0.07. 



A. 38. GRB 080805 

The optical spectrum of the afterglow of GRB 08080 5 (z = 
1.5042 ) was obtained with the VLT/FORS2 (Jakobsso rietall 
'2OO84 F09). The spectrum clearly show a 2175 A absorption 
feature at the redshift of the GRB. Due to the relative faintness 
in the optical, the afterglow is not detected in any of the UVOT 
filters. The spectrum shows a pretty clear absorption dip and the 
extinction curve resembles that of the MW. The intrinsic SED of 
the afterglow is generated 1.0 hr after the burst trigger The X- 
ray to optical SED is reproduced nicely with a single power-law 
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plus FM extinction curve with a large amount of dust redden- 
ing of Ay = 1.53 + 0.13 mag and total-to-selective extinction of 

Rv = 2.48 + 0.39^ 

Previously iGreiner et al ] (12011b modeled the GROND-XRT 
data and found that MW dust model, with a 2175 A bump, pro- 
vides the best explanation to the data with Ay - 1.01^y,'g. The 
value is consistent with the results we find from our analysis. 

A.39. GRB080905B 

The optical spectroscopy of the afterglow of GRB 080905B 
at z = 2.3739 was taken with the VLT/ FORS2 together 
with th e R band acquisition image photometry dVreeswiik et al.l 
l2008bt F09V We implem ented MIDAS package FITLYMAN 
jFontana & Balles"te^ll995h to determine the neutral hydrogen 
column density. The derived H i column density of the afterglow 
of GRB 080905B is log Nui/cm-^ < 22.15. Since the afterglow 
flux is contaminated by the light from another object on the slit, 
we suggest the column density to be an upper limit. We modeled 
the intrinsic afterglow SED and the data are well fitted with a 
single power-law and dust extinction with Ay - 0.41 + 0.03 mag 

A.40. GRB 08091 3 

The spectrum of t he afterglow of GRB 0809 13 was taken with 
the VLT/FORS2 dGreiner et al.ll2009t F09: IZafar et"an 1201 lb . 
GRB 080913 is the highest redshift GRB in our spectroscopic 
sample. The acquisition image was taken with z-Gunn filter. The 
magnitude reported in Table [1] is strongly affected by Lya ab- 
sorption. We corrected the magnit ude for Lya absor ption using 
spectro-photometric analysis (see IZafar et al.l 1201 ll for detail). 
A comprehensive description of the SED analysis is given in 
IZafar et al.l (1201 ih . T he neutral H i colurn n density of tiie DLA 
is Nm/cm-^ < 21.14 (iGreiner et al.ll2009l) . 

A.41. GRB080916A 

The optical spectrum of the optical afterglow of G RB 0809 16A 
(z ^ 6887) was carried out with the VLT/FORS2 dFvnboetalJ 
l2008al : F09). The SED is constructed at 17.1 hr i.e. the time ofR- 
band acquisition image photometry. The intrinsic SED is nicely 
reproduced with a broken power-law and SMC type dust extinc- 
tion with Av - 0.15 + 0.04 mag. 

A.42. GRB 080928 

The optical spectrum of the afterglow of GR B 080928 at z - 
1.6919 was obtained with the VLT/FORS2 dVreeswiik et al.l 
l2008at F09). We collected the optical/NIR photometry at fli e 
time of the SED taken with GROND from Rossi et al. (1201 ll) . 
We scaled the spectrum to the level of the r' band photome- 
try taken at 15.5 hr. In addition we reduced the UVOT data for 
this burst and fluxes are obtained for the v, b and u bands at 
the common SED epoch. The intrinsic SED is well fitted with a 
single power-law and SMC extinction. The best fit Ay value is 

0.3 +0.03 mag . 

iKann et alj d2010l) find that a MW dust model best fits the 
optical-NIR data with a small amount of extinction (Ay - 0-14 + 
0.08). We do not find any evidence for the 2175 A bump in the 
optical spectrum. 
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Fig.A.l. Afterglow SEDs of the spectroscopic GRB sample in vFy and v space. In each figure, we plot the best fit extinguished 
(solid lines) and extinction corrected spectral model (dashed lines). The grey curve represents the optical spectrum. The X-ray 
spectra are indicated by blue points. The black triangle corresponds to acquisition camera photometry used for scaling the SED. 
Black open circles are not included in the spectral fitting because of the optical spectrum wavelength coverage in that region while 
solid circles represent the data points included in the SED modeling. The red curve shows the best fit single power-law model and 
the green curve illustrates the broken power-law model. 
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Fig. A.l. (continued) 
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Fig. A.l. (continued) 
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Fig. A.l. (continued) 
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Fig. A.2. Absolute extinction curves of the spectroscopic GRB sample. The extinction curves are based on the best fit models given 
in Table |2l The grey curve represents the optical spectrum. The black triangles correspond to acquisition camera photometry used 
for scaling the afterglow SED. Black open circles are not included in the spectral fitting because of the optical spectrum wavelength 
coverage in that region while solid circles represents the photometric data points included in the SED modeling. The red solid curve 
corresponds to the best dust mod e l for e ach GRB. Also shown are the Milky Way (green dashed line) and SMC (blue dotted line) 
extinction models taken from |Peil(ll992h . 
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Fig. A.2. (continued) 
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